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ABSTRACT 



Context, X-rays from very young stars are powerful probes to investigate the mechanisms at work in the very first stages of the star 
formation and the origin of X-ray emission in very young stars. 

Aims. We present results from a 500 ks long observation of the Rho Ophiuchi cloud with a XMM-Newton large program named 
DROXO, aiming at studying the X-ray emission of deeply embedded young stellar objects (YSOs). 

Methods. The data acquired during the DROXO program were reduced with SAS software, and filtered in time and energy to improve 
the signal to noise of detected sources; light curves and spectra were obtained. 

Results. We detected 111 sources, 61 of them associated with p Ophiuchi YSOs as identified from infrared observations with 
ISOCAM. Specifically, we detected 9 out of 11 ClassI objects, 31 out of 48 Class II and 15 out 16 Class III objects. Six objects 
out of 21 classified Class III candidates are also detected. At the same time we suggest that 15 Class III candidates that remain unde- 
tected at logL x [erg/s] < 28.3 are not related to the cloud population. The global detection rate is ~ 64%. We have achieved a flux 

sensitivity of ~ 5 • 10~ 15 erg s cnT 2 . The Lx to Lboi ratio shows saturation at a value of 3.5 for stars with T e u < 5000 K or 0.7 Mo 

as observed in the Orion Nebula. The plasma temperatures and the spectrum absorption show a decline with YSO class, with Class 
I YSOs being hotter and more absorbed than Class II and III YSOs. In one star (GY 266) with infrared counterpart in 2MASS and 
Spitzer catalogs we have detected a soft excess in the X-ray spectrum, which is best fitted by a cold thermal component less absorbed 
than the main thermal component of the plasma. This soft component hints at plasma heated by shocks due to jets outside the dense 
circumstellar material. 

Key words. Stars: activity. Stars: formation. X-rays: stars. X-rays: individuals: Rho Ophiuchi 



1. Introduction 

Low-mass stars in the pre-Main Sequence (PMS) phase are char- 
acterized by intense X-ray emission. X-ray activity has not been 
yet firmly confirmed in the very initial protostellar cores that 
accrete material from a surrounding envelope (Class objects), 
because X-rays are either absent or completely absorbed by the 
dense circumstellar material. In Class I, II, and III objects and 
until the Zero Age Main Sequence stage the X-ray emission is 
very strong when compared to solar-age stars and also charac - 
terized by impulsive variability dFeigelson & MontmerleL [l999). 
In the case of Class III objects all X-rays are thought to origi- 
nate from a magnetized stellar corona, which likely resembles 
a scaled-up version of the corona of late-type Main Sequence 
stars and the Sun. Additional X-ray production mechanisms are 
possibly at work in Class I and II objects. At these young evo- 



lutionary stages the magnetic field drives the material in-falling 
from the disk to the stellar surface where the matter becomes 
shocked. Furthermore, outflows along the disk axis may interact 
with the circumstellar material and form shocks. The plasma in 
the shocks in both in- and outflows may be come heated up to a 
few million degrees, thus emitting X-ray s (Pravdo et al., 2001; 
iFavata et all l2002t iKastner et all 120051 iGiardino et ail [2007; 



Isacco et all l2008). 
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Significant effort has been devoted to the X-ray observa- 
tion of star-forming regions, which provide natural laborato- 
ries where to study the X-ray emission from young stars and 
its implications for the mechanisms of star formation. The p 
Ophi uchi cloud is amon g the nearest star-forming regions (120 
pc. lLoinard et all 12008) and has been extensively studied from 
the infrared (IR) to the X-rays bands. The p Ophiuchi cloud is 
shaped as a dense multicore structure, hosting more than 200 
members comprising young stellar objects (YSOs) in all evo- 
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lutionary stages from Class to Class III ( Wilki ng et"aTl 12008). 
While a dispersed population of optically visible stars associated 
with the Upper Scorpius OB association with an age of ~ 5 Myr 
is present in the region, the studies in the IR band have revealed a 
younger population of Pre Main Sequence (PMS ) stars and pro- 
tostel lar objects with ages of only 0.3-1 Myr dLuhman & RiekeL 
1 19991) . which are thus younger than YSOs in other star-forming 
regions like the Taurus Mol ecular Cloud (~ 1 - 5 Myr) and 
the Orion Nebula (~ 1 Myr, iHillenbrandl 1 19971) . The mid-IR 
survey with ISOCAM on board the ISO satellite reported by 
Bontemps et al.l d2001l) (hereafter BoO 1 ) allowed the detection of 
a population of 16 Class I, 123 Class II and 77 Class III YSOs, 
adding 71 previously unknown members of Classes I and II. X- 
ray studies carried out with the Einstein and ROSAT satellites 
had revealed several tens of embedded Class I, II, and III stars 
that are highly active i n X-rays and are charac t erized by a very 
strong time variabil ity ( Montmerl e et all 1 1 9831; ICasanova et al.L 
119951) ICasanova et al.l (1 1995b tentatively identified several X-ray 
sources with Class I protostars with ROSAT/PSPC observations. 
However, X-ray observations with the ROSAT/High-Resolution 
Imager we re needed to confirm X-ray emi ssion from Class I 
protostars dGrosso et all 119971: iGrossoL 1200 1|). Re cent observa- 
tions with the XMM-Newton dOzawa et all 2005) and Chandra 
satellites dlmanishi et all l2001[|Flaccomio et al mool have in- 
creased the number of candidate PMS members of the region. 
These studies suggested that accreting YSOs have increasing X- 
ray activity going from Class I to Class III and decreasing plasma 
temperatures. The aim of the Deep Rho Ophiuchi XMM-Newton 
Observation (DROXO) was to obtain a high-sensitivity survey 
in the Core F of p Ophiuchi Cloud by means of a long, almost 
continuous observation. We report the data analysis and the X- 
ray properties of the PMS stars observed during the DROXO 
program. The structure of the paper is as follows: Sects. 2 and 3 
describe the observation and the data analysis, Sects. 4 and 5 de- 
scribe the sensitivity of the survey and the nature of the sources 
detected in DROXO. In Sect. 6 we discuss the X-ray emission of 
the sample of classified YSOs, in particular we explore the be- 
havior of plasma temperatures, absorption, X-ray luminosity and 
ratio Lx/Lboi with respect to effective temperatures, masses and 
evolutionary stages. In Sect. 7 we give a summary. Appendix 
A and B contain tables with the list of detected sources, the 
best-fit parameters of spectra, the upper limits to rate for un- 
detected YSOs in the ISOCAM sample, the Spitzer counterparts 
to DROXO sources, and an example of one page of the atlas 
(published electronically only) with the EPIC spectrum and light 
curve of each source. 



2. The observation 

The Deep Rho Ophiuchi XMM Observation is a large program 
carried with the XMM-Newton satellite pointed toward core 
F in the p Ophiuchi Cloud for almost eight consecutive days. 
The nominal pointing was at R.A. = 16h27ml9.5s and Dec. = 
-24d41m40.9s (J2000), and the net exposure time was ~515 
ks. The pointing approximately co incides with that of the 33 ks 
XMM-Newton exposure studied bv lOzawa et al.l d2005h . The ob- 
servation has been carried out in five subsequent satellite orbits 
(0961 — » 0965), keeping the same position angle in all orbits. 

During the first orbit chip nr. 6 of the MOS 1 camera 
was damaged apparently by a micrometeorite impact Q and has 
stopped functioning, so that only a ~ 28.5 ks exposure is avail- 
able by that chip in DROXO. No damage was registered on the 



other chips of MOS 1 . During orbit nr. 0964 the p Ophiuchi ex- 
posure was split in two segments, separated by ~ 25ks, due to a 
triggered TOO observation. During the first three orbits the drift 
of the satellite was larger than 6", and this influenced the compu- 
tation of exposure maps and the subsequent analysis of summed 
data as explained below. 

In Fig. Q] we show the EPIC MOS1, MOS2 and PN images 
added in a pseudo-color Red-Green-Blue image. The bands cho- 
sen for red, green and blue components are: 0.3-1.0 keV, 1.0-2.5 
keV, 2.5-8.0 keV. Sources with a softer/less absorbed spectrum 
are redder than sources with a hard/highly absorbed spectrum. 
Each CCD image has been divided by the proper exposure map 
and by the average effective area in the energy band to normalize 
the different efficiencies of the three cameras. 



3. Data analysis 

The observation data files (ODF) were processed with the SAS 
softwar^l (version 6.5) to produce full field-of-view event lists 
calibrated in both energy and position. We subsequently filtered 
these event files and retained only the events with energy in 
the 0.3-10 keV band and those that triggered simultaneously at 
most two nearby pixels. This step was executed for each detec- 
tor (MOS 1, MOS 2, PN) and for each of the exposure segments 
of the five orbits (hereafter step 1). The large satellite drift is 
not taken into account automatically by the SAS software which 
results in wrong exposure maps when they are created with the 
default values. We needed to reduce the ATTREBIN parameter 
of task EEXPMAP to 0.5 " and to choose the most accurate al- 
gorithm (of which the default was the fastest and less precise 
one) to obtain correctly computed exposure maps. 

3.1. Source detection 

We performed source detection with the PWXDETECT code 
developed at I NAF-O sservatorio Astronomico di Palermo 
' fl997bllah . The code allows the detection of 



(Da miani et al.l 

sources starting from unbinned photon positions recorded in sev- 
eral datasets from MOS 1, MOS 2 and PN cameras, through a 
multiscale analysis of mexican-hat wavelet convolved images. 
All data were properly scaled by time and effective area of each 
CCD detector, obtaining a flux image of the sum of all expo- 
sures. At the end of the process the count rates of detected 
sources were re-scaled to a reference instrument which, in our 
case, was the EPIC MOS 1 . We refer to the count rates derived 
in this way as MOS1 equivalent count rates. In order to improve 
the detection of faint sources, we filtered the data obtained at step 
1, excluding those time intervals with a high background count 
rate. In fact, the background during DROXO was highly variable 
for a significant fraction of the exposure time (see Fig. [2]). We 
excluded all events registered in intervals when the rate on the 
whole image, shown in Fig. [2] was higher than a given thresh- 
old. The rate threshold was chosen in a way that maximizes the 
signa l -to-noise-ratio (SNR) for faint sources (cf. Sci ortino et all 
2001; Damian i et al 1 120031) and improves the source detection 
process toward faint sources as expected. The net cleaned expo- 
sure times after this filtering are 198.1 ks (38%) 273.2 ks (53%) 
and 213.2 ks (41%) for MOS 1, 2 and PN, respectively. 

We detected 111 point sources with a significance thresh- 
old corresponding to two spurious detections in the whole field- 
of-view. The threshold for the detection of faint sources was 
determined from the analysis of a large set of simulations of 



1 http://xmm. vilspa. esa. es/external/xmmjiews/items/MOSl - CCD6/index. shtinbee http://xmm. vilspa. esa. es/sas 
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Red: 0.3-1.0 keV; green:1 .0-2.5 keV; blue:2. 5-8.0 keV 
H ' ' ' ' I I ' ' ' ' H 




+ 16:28:0.0 +16:27:30.0 +16:27:0.0 +16:26:30.0 

R.A. (J2000) 

Fig. 1. Merged EPIC image of events recorded in the time- 
filtered data to enhance the signal-to-noise of faint sources. 
Colors encode the following bands: 0.3-1.0 keV (red), 1.0- 
2.5 keV (green), 2.5-8.0 keV (blue), respectively. The MOS 1, 
2, and PN images are normalized by effective area and exposure 
time to reduce instrumental artifacts like CCD gaps. 



background-only images and then running the detection code on 
these images. The simulations provide a value for the signifi- 
cance threshold to retain at most 1-2 spurious sources per field 
in real data. 

Source positions, off-axis distance, exposure times, and X- 
ray count rates are listed in Table Al for all detected sources. We 
also report the 2 MASS designation and other names from the 
literature for the optical counterparts identified in Sect. 5.1. The 
last column of Table Al contains a flag pointing to the source of 
the previous identifications. We built an atlas of spectra and light 
curves for each source which is available onlinqj; in Appendix 
|B]we show a page of this atlas as an example. 

3.2. X-ray spectra 

We produced both light curves and spectra of all sources by se- 
lecting the photons in circular regions around the source posi- 
tions. The radii of the regions depend on the source intensity, 
on possible of nearby sources and the geometry of the CCD. 
We used regions on the same CCD for source and background; 
for the background, we avoided to include pixels and columns 
poorly calibrated in energy. Figure [3] shows two examples of 
the choice of source and background extraction regions. For PN 
spectra a further constraint was to have both source and back- 
ground extraction regions at nearly the same distance from the 
CCD readout node. 

To produce the spectra we filtered the photons with respect to 
FLAG (chosen to be equal to zero) and PATTERN (less or equal 
to four) as indicated in the SAS User's Guide. Given that the 



3 See: http:www.astropa.unipa.it/ pilli/ atlas_droxo_sources.pdf 
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Fig. 2. Top panel: light curve of all events recorded with MOS 
1. Labels for the individual satellite orbits are given on top of 
the plot. The gray shaded area represents the time interval when 
MOS 1 was turned off after the micro-meteorite impact. The hor- 
izontal line is the threshold count rate that maximizes the SNR 
of the full image as described in the text. Bottom panel: signal- 
to-noise function vs cumulative time. 



choice of GTIs was tailored toward the faintest sources, at this 
stage we further improved the choice of GTIs for bright sources. 
The time-filtering we adopted for spectra starts from the Good 
Time Intervals (GTIs) defined initially to perform source detec- 
tion and iteratively adds temporal bins from the light-curve Q 
that contribute to increase the total SNR of the spectrum. The 
procedure starts by adding the time bin with the largest individ- 
ual SNR and is then iterated considering time-bins of decreasing 
SNR until no gain in total SNR is obtained. As a result, the net 
exposure time of each spectrum is different. 

3.2.1. Background subtraction 

We noticed that scaling the background by the geometric areas of 
the source and background extraction regions leads to incorrect 
estimates of the background, especially for faint sources and/or 
during times of high background. The effect is that background- 
corrected light curves of faint sources are either directly or in- 
versely correlated with the background light curves and often 
lead to negative net count-rates. We understood this as the ef- 
fect of a vignetted background component that is not properly 
taken into account by a purely geometric scaling factor. In order 



4 The binning of the light curve for each source is chosen in an itera- 
tive way, starting from Ar = 1 ks and progressively enlarging AT until 
the bin with maximum net counts has at least 50 counts or AT becomes 
larger than 10 Ks. 
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Fig. 3. Two examples of choices for the source and background extraction regions. The left panel shows two nearby but distinct 
sources. The right panel shows an extreme case in which the region for the faint source is shrunk to minimize the influence of the 
bright source. The background regions are chosen as near as possible to the sources to which they refer and at the same distance of 
the readout node for PN. 



to correctly estimate the background contribution to the photons 
extracted in the source region we constructed background maps 
for each instrument (MOS 1, 2 and PN) and for each orbit. These 
were built by removing large regions around detected sources 
from the images, and by subsequently smoothing and interpo- 
lating the maps over the source extraction areas. Although the 
resulting scaling factors differ from the purely geometric ones 
by only a few percent, the difference is relevant in cases when 
the background dominates the count rate in the source regions 
and mitigates the above mentioned spurious effects on the light 
curves. 



3.2.2. Model fitting of spectra 

We analyzed the spectra of all X-ray sources in the 0.3-10 keV 
band. For a given EPIC detector (MOS 1, MOS 2 or PN) spec- 
tra from all orbits were summed up; analogously, background 
spectra were obtained; the response matrices and ancillary re- 
sponse files for each spectrum were multiplied and then summed 
by weighting by the exposure time. The background was scaled 
according to the procedure described in Sect. 13.2.11 

The spectra were grouped prior to fitting with XSPEC v. 12.3 
to obtain at least a minimum SNR in each bin, by consider- 
ing both source and background photons. In order to obtain the 
largest number of meaningful grouped spectra we adopted two 
schemes of grouping procedure based on high and low SNR of 
the final spectrum. For this purpose we used the same procedure 
as in the ACIS_EXTRACT package for the analysis of Chandra 
ACIS datcQ, adapted to our EPIC spectra to take into account the 
background. 

We grouped the spectra using two thresholds for the mini- 
mum SNR to be obtained in each bin, i.e. we generated two sets 



s See http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
for the details of the algorithm. 



of spectra, one set with a high and one with a low SNR per bin. 
The minimum SNR in each spectral bin was imposed on the ba- 
sis of the source count statistics, SNR > 2 defining our low SNR 
binning and SNR ~ 3 - 5 defining our high SNR binning. Where 
possible, we tried to obtain binned spectra with at least eight 
bins. 

The spectra of all detectors were fitted simultaneously. In 
some cases we had to discard the data of one or two of the three 
EPIC cameras. These cases occur when a source is on CCDs 
gaps of one or two cameras, thus leading to a wrong estimate 
of the point spread function fraction contained in the extraction 
area. 

The spectra were fitted with one-temperature (1-T), two- 
temperature (2-T) and, in some cas es, three temperature (3- 
T) APEC models (ISmifh et all 1200 ll) plus absorption (WABS) 
dMorrison & McCammonL 1 1983h . the free parameters were the 
absorption column Nh, the temperatures, and the emission 
measures. The abundance pattern was fixed to that found 
in P MS coronae f rom t he Orion Nebula in the COUP sur- 
vey dMaggio et all 120071) . Only in four cases described below 
(Elias 29, SR 12, IRS 42, src. 61) a more complex model was 
required. 

For Elias 29 (src. 38) and SR 12A (src 53) the global abun- 
dance scaling was left as a free parameter to achieve a bet- 
ter fit. For Elias 29 we confirm an unusually high abundance 
(Z ~ Z ) with respect to those found in PMS co ronae (Z ~ 
0.2 - 0.3Z Q ), as already reported by iFavata et all (120051) . The 
source IRS42/GY252 (src. 54) is located in the wings of the 
much brighter X-ray source corresponding to SR 12A (see right 
panel of Fig. 3). To take account of the contribution by SR 12A 
we added the best fitting 3-T APEC model multiplied by a con- 
stant factor representing the amount of contamination to the 
model of IRS 42/GY252. This scaling factor was a free fitting pa- 
rameter. The spectrum of IRS 42/GY252 itself can be described 
with a 1-T APEC model. The spectrum and the best-fit model 
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Fig. 4. Spectrum and best-fit model of IRS42/GY252. The 
stepped curves are the contributions of SR 12A (light line) and 
that of IRS42 (dark line). The spectrum of IRS42/GY252 is 
modeled by an absorbed single temperature APEC model. 



are shown in Fig. [4] As can be seen, IRS42/GY252 is highly 
absorbed and the soft emission is attributed entirely to the con- 
tamination by SR 12 A. For the source nr. 61 we used a model 
given by the sum of two differently absorbed APEC models to 
account for the soft excess visible below 1 .0 keV, as discussed 
in Sect. I6\4land Fig. [TT1 

Generally, the best fit was chosen on basis of the probabil- 
ity P(x 2 > Xq) of obtaining a x 2 higher than the observed one 
Xq- Our threshold for acceptable fits was P(x 2 > xfy > 0.02. 
Whenever this criterion was satisfied by the spectra with high 
SNR per bin (typically for medium-high statistics spectra) we 
chose these; otherwise we selected the fits obtained from the 
spectra with low SNR per bin. We always selected the best-fit 
model with the smallest number of free parameters. In seven 
cases, even for some bright sources, no formally acceptable fit 
was found and we allowed for a lower P(x 2 > Xq)- However, we 
accepted the best fit results by visually checking that the overall 
shape of the spectrum is well reproduced by the model. 

With this procedure we obtained spectral fits for 91/111 
sources. Table [A2l summarizes the results. The columns report 
the source number, the available EPIC datasets, the data sets we 
chose for fitting, a flag indicating the binning of the spectra we 
used (high or low SNR), the type of model, the model best-fit pa- 
rameters, the unabsorbed flux and luminosity in the 0.3-10 keV 
band, the x 2 statistics, degrees of freedom and the probability 
P(X 2 > Xq)- F° r those weak sources without a spectral analysis, 
we calculated fluxes and X-ray luminosities in the 0.3-10.0 keV 
band by using their count rates and PIMMS softwar^] assuming 
a 1-T model with kT and N# equal to the median of values de- 
rived from the best-fit procedure of spectra (see Sect. HJ and a 
distance of 120 pc. 

3.3. Photospheric parameters 

In Sects. 5 and 6 we will focus our discussion on the X-ray prop- 
erties of the sample of PMS stars classified by BoOl. The photo- 
spheric stellar parameters (r e ff, L\, \ and masses) for the Class II 
and III objects from the list of BoOl were estimated from the 



6 see http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



near IR (2MASS) photometry. The procedure we u sed closely 
follow s that adopted by BoOl and improved by iNatta et akl 
(2006). We assumed that the J-band emission from these sources 
is dominated by the stellar photosphere and is only marginally 
contaminated by the emission from circumstellar material and 
also that the IR colors of Class II sources can be described as 
the emission from a passive circumstellar disk as described by 
iMever etaTI £l997). These assumptions obviously do not apply 
to Class I sources and for this reason photospheric parameters 
for these objects were not estimated. 

We used the ICardelli et ail (Il989l) extinction law with R,, = 
4.4, which we think is appropriate for Ophiuchus. A small num- 
ber of sources (~15%) have colors slightly bluer than those of 
reddened main sequence stars, presumably due to photometric 
uncertainties, which are on the order of 0. 1 mag, while the off- 
sets of these objects with respect to the reddened sequence range 
between and 0.15 mag. Dereddening these sources by extrap- 
olating the colors of Class II and III sources would produce an 
overestimate of the extinction. For these objects we assigned the 
colors of the closest photosphere model on the reddened main 
sequence. Table lA.5l lists the effective temperatures, masses and 
bolometric luminosities of ISOCAM objects. 

The values of t he 7-band extincti on we derive are very sim- 
ilar to ones from INatta et a"D d2006l) . with only the significant 
exception of WL 16, for which our procedure produces a signif- 
icantly higher extinction. 



4. Sensitivity of the survey 

The faintest detected source in DROXO has a net count rate of 
8.3 • 10~ 2 ct ks~' in MOS equivalent units. In Fig. [5] we plot 
a sensitivity map in units of count rate per point source of the 
DROXO field of view. It is built starting from a smoothed back- 
ground map, after removing the contributions of the sources, and 
taking into account the vignetted exposure map and the thresh- 
old used for detection. The map shows that the sensitivity varies 
by a factor 2.5 in the area covered by the three EPIC cameras, 
and it is quite constant across a 60% fraction of the field of view 
because of the characteristics of the detector. To translate the 
limit count rate in a limit flux we needed a conversion factor 
that depends on the spectrum temperature and, critically, on the 
absorption. From the analysis of the spectra we know that the ab- 
sorption is also strongly variable from source to source by more 
than a factor 100, which can be due to local material around the 
sources or dense cloud material in front of the objects, or both. 
The plasma temperatures that we find from the spectra vary from 
kT 55 1 keV to a few ke V. The combined action of local absorp- 
tion and plasma temperatures of undetected sources causes the 
limit flux to vary more than the range observed in the sensitivity 
map showed in Fig. [5] 

For example, we can estimate a conversion factor from count 
rates to fluxes by using the median of the column absorption 
(Nh = 2.3 • 10 22 cm~ 2 ) and of the plasma temperatures (kT = 3.1 
keV) from the spectral fits, and this gives us cf = 6.1 • 10~ u erg 
cts~' crrT 2 . This yields a limiting unabsorbed flux in the 0.3- 
10 keV band of ~ 5 ■ 10~ 15 erg s _1 cirT 2 and a luminosity of 
~ 8.7 ■ 10 27 erg s _1 . The c/s derived from all sources are in the 
range cf = 2.3 - 12.7 ■ 10~ n erg cm~ 2 cts -1 , yielding a limit 
flux comprised in the range Fx = 5- 22 -10-15 erg cirT 2 s _1 
and luminosities Lx — 7 - 38 • 10 27 erg s , respectively. These 
values are indicative of the sensitivity we achieved across the 
field of view, but local strong absorption can significantly lower 
the actual sensitivity in that position. 
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Fig. 6. Left panel: MOS 1 rate from lOzawa et al. I (12001 vs. DR OXO MPS 1 rate . The lines mark the region for rate variations 
of factors 2 and 5. Right panel: the same for the ACIS rate from Flaccomi o et al.l d2003l) . Upper limits to count rate of sources 
undetected in the compared samples are indicated by vertical and horizontal arrows. Errors are quoted at l<x level. 
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Fig. 5. Sensitivity map of the DROXO field in units of ct / s per 
point source (logarithmic scale). The map refers to the sum of 
PN and MOS instruments. The intensity scale is indicated in the 
vertical strip on the right side of the map. Crosses indicate the 
positions of detected sources. 



4.1. Comparison with other X-ray surveys of p Ophiuchi 

In Fig. |6l we s how t he scatter plots of the MOS 1 count rate from 
lOzawa et al.l (120051) and Chandra- ACIS rate (wh ich has an effec- 
tive area comparable to the EPIC MOS) from iFlaccomio et al.l 



(120031) vs. the DROXO MOS equivalent rat e. In DROXO we de - 
tected 80 of the X-ray sources reported by Oza wa et al.l ([2005 ), 
while 7 of the sources found by Ozawa and collaborators remain 
unde tected in DROXO; ana logously we identified 62 sources 
from [Fl ac comioetaL (2003), while 30 of those sources have not 
been detected. We computed an upper limit to the cou n t rates 
for those sources that a re detected by IFlaccomio et al. I J200l 
and lOzawa et all d2005) but are undetected in DROXO (indi- 
cated as arrows in Fig. [6). We also re-analyzed the pre vious 
XMM-Newton observation reported bv lOzawa et al. (2005) with 
the same detection procedure as in DROXO, and we calculated 
the upper limits to count rates for DROXO sources undetected 
in the Ozawa survey. Analogously we calculated upper limits to 
co unt rates for two sourc es detected in DROXO but undetected 
bv IFlaccomio et all (120031) . The count rates of sources detected 
in DROXO and in the two other surveys globally agree within a 
factor ~5. Because there is no systematic trend in the scatter we 
conclude that it is likely to be attributed to variability. Although 
a detailed study of the time variability is beyond the scope of this 
paper, we report here that 52% of the sources have variable rates 
at more than 90% significance level when compared with the 
Ozawa et al. survey, and 79% are found variable to be compared 
with the survey of Flaccomio et al. 

In the 33 ks XMM-Newton survey of lOzawa et all d2005l) the 
faintest detected source has a count rate of 0.5 ct ks~'. Using 
the conversion factor derived by us for DROXO the correspond- 
ing luminosity is 5. 3 TO 28 erg s _1 , which scales reasonably well 
with the exposure times of both surveys. Flaccomio et al.'s sur- 
vey reaches a limiting rate similar to DROXO, likely due to the 
lower Chandra-ACIS background with respect to that found in 
our EPIC observation. Most of the ACIS sources undetected in 
DROXO are near the sensitivity limit of DROXO, and a variabil- 
ity of factor 2 can easily explain their missed detections. 

Sources 33 (GY 304) and 43 are brighter in DROXO than in 
Ozawa et al. survey by more than a factor five. DoAr 25 suffered 
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Time (ks) after 2005-03- 1 6T22:1 



Fig. 7. PN Light curve of WL2 (source nr. 15). The bold light 
curve is the net rate, while the thin line is the total rate (source 
+ background). A large flare lasting ~35 ks is well visible. The 
source quiescent rate is comparable to the background rate. 



of strong pile-up in the Ozawa et al. survey. These authors do 
not report a direct measurement of the count rate for this object, 
but they derived a luminosity from the spectral analysis. From its 
X-ray luminosity we inferred a count rate greater than 200 ct/ks. 
We checked that in DROXO DoAr 25 and the other two brightest 
sources in the field, SR 12A and IRS 55, are not affected by 
pile-up. The variability of DoAr 25 with respect to the XMM 
observation detailed in Oza wa et a l] (l2005h is at least a factor 4. 

Among the most variable sources, WL 2 (labeled in Fig. [6]l 
shows a big flare in the DROXO light curve (Fig. |7). The qui- 
escent rate of WL2 is ~ 6 ct ks~' and the peak rate is ~ 300 
ct ks~' explaining its large offset position in the scatter plot of 
the right panel in Fig|6] The quiescen t rate is consistent with the 
one measured bv lOzawa et al.l d2005l) . The long duration of this 
flare on WL2 (~ 35 ks) underlines the need for long observa- 
tions of YSOs to properly assess their quiescent emission. 



5. The nature of the X-ray sources 

5. 1 . Optical and IR counterparts of the DROXO sources 

We searched for optical and IR counterparts of D ROXO sources 
in the 2MASS, Spitzer C2D dEvans et alll2003alfl and ISOCAM 
(BoOl) surveys a s well as in the opti c al and IR surveys of 
iNatte et all (f2006l). IWilking et all (l2005h . iBarsonv etail (120051) 
and iLuhman & Rieke (Il999t) . The match radius takes into ac- 
count the uncertainties on the X-ray source position in DROXO, 
which is generally on the order of ~ 1" , and other catalogs. The 
mean value for the match radius is ~ 1.0" with a 0.1-0.9 quan- 

tile range of ~ 0.5 5.8". As anticipated in Sect. 3.1, Table 

I A. II reports the literature names of DROXO counterparts and a 
flag indicating in which catalogs the source is identified. Sources 



with counterparts in the ISOCAM survey are indicated with let- 
ter "I" in the "Flag" column; sources identified only in Spitzer 
are indicated with letter "S" ; s ources which hav e only a X-ray 
identification in the surveys by [Flaccomio et al.] d2006l) and/or 
Oza wa et a D d2005h are indicated with letter "X". Among the 
111 X-ray sources 78 have a counterpart in 2 MASS and 61 have 
a counterpart in ISOCAM data. Two objects have no 2 MASS 
and no ISOCAM counterpart, but have been identified in the lit- 
erature (IRS50 and WSB46). 

An additional 16 X-ray sources have been newly identified 
with Spitzer. This leaves 16 DROXO sources without known op- 
tical/IR counterpart. Of these, 6 were detected in previous X-ray 
surveys and 10 are presented here for the first time. We calcu- 
lated upper limits to the X-ray fluxes and luminosities for the 
sample of undetected ISOCAM objects falling in the field of 
view of DROXO using the conversion factor derived in Sect. [4] 



5.2. Unidentified X-ray sources 

The inspection of light curves and spectra can give some indi- 
cation on the nature of sources with unknown optical/IR coun- 
terpart, designated 'U' or 'X' in Table I A. 1 1 Five unidentified 
sources (Src. # 5, 7, 16, 26, 58) are too faint for spectral anal- 
ysis, but the other five unidentified X-ray sources (Src. #12, 
19, 88, 96, 110) have spectral parameters compatible with the 
expectation for a YSO in pOph, i.e. kT of a few keV and 
Nn > 10 22 cnT 2 . Furthermore, Src. 110 shows impulsive time 
variability similar to flares typical for PMS stars. There are also 
six sources (Src. # 20, 29, 41, 45, 48, 95) without optical/IR 
counterpart that have been detected in previous X-ray surveys. 
For Src. # 20, 29 and 41 we had too few counts to obtain mean- 
ingful spectral fits, and their light curves show some variability, 
but no clearly identifiable flare. For Src. # 45, 48 and 95 the spec- 



tral analysis gives a high absorption (above 10 22 cm^). While 
the temperature of Src. # 95 is not constrained, Src. # 45 and 
48 have plasma temperatures of 5.4 and 4.5 keV. For these two 
sources a power law best fit to their spectra is also acceptable. 
These three objects could have characteristics consistent with 
those of highly embedded YSOs, although they are undetected 



See also http://ssc.spitzer.caltech.edu/legacy/c2dhistory.html 



on the millimetric survey s of Motte et al.ldl998l) : lJohnstone et al.l 
(120001) ; |j0rgensen etail d2008l) . The lack of 2MASS counter- 
parts suggests that these sources could have an extragalactic na- 
ture. On the other hand, it is not ruled out tthat they are very low 
mass PMS stars or even brown dwarfs. 



5.3. Mid-IR photometry from Spitzer 

The Spitzer photometry for the X-ray sources identified in the 
C2D survey are listed in Table A4. We derived magnitudes in 
the 3.6, 4.5, 5.8, 8.0 /.im IRAC bands and constructed the color- 
color diagram shown in Fig.[8](color index [3.6] - [4.5] vs color 
index [5.8] - [8.0]). This diag ram provides a rou gh classification 
of YS Os according to Alle n et al. I d200l and lHartmann et al.l 
d2005h . Normal, unreddened stars or Class III / Weak T-Tauri 
stars with very low reddening should occupy a region centered 
around (0,0). Reddening due to matter along the line of sight 
tends to disperse vertically data points along the [3.6]-[4.5] 
IRAC color index, while the [5.8]-[8.0] color index is not af - 
fected by this source of reddening (see Flahe rty et all 120071) . 
Infrared emission from a circumstellar disk in Class I and II 
sources moves the objects both vertically and horizontally to- 
ward redder values of both indexes. We expect to find Class 
II YSOs in the region marked with the black box or above it 
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Fig. 8. IRAC color-color diagrams of objects in the DROXO field with information on the reddening from J or V band photometry 
or from A^h- The left panel shows colors before correction for the reddening, the right panel shows dereddened IRAC colors. Purple 
squares mark objects with X-ray emission in DROXO. Asterisks are Class III YSOs, triangles are Class II YSOs and diamonds are 
Class I YSOs. The rectangle marks the region where Class II YSO colors are expected. 



and very embedded protostars (Class 0/1) are expected to lie at 
[3.6] - [4.5] > 1.1 and at [5.8] - [8.0] > 1 (lAllen et al.ll200l . 

We plotted only the objects for which we have information 
on the reddening from J and V band photometry or from A^h- 
In the left panel of Fig.[8]the symbols represents the IRAC col- 
ors before correction for reddening. Objects detected in DROXO 
are marked with squares, diamonds are Class I YSOs as classi- 
fied by BoOl, triangles are Class II YSOs, asterisks are Class III 
YSOs. The points with ISOCAM classification follow the distri- 
bution outlined above, with Class III YSOs around or dispersed 
vertically above the (0,0) point, and Class II YSOs located in the 
rectangle or dispersed vertically above this region. 

In order to investigate the nature of these latter "dispersed" 
objects we dereddened the colors of those objects for which we 
had information on the extinction. If available we used the red- 
dening Aj or Ay from the literature, otherwise the column ab- 
sorption jVh from the X-ra y spectra with the conversion factor 
of A v = 1.51 • 10~ 22 ■ N A dMathisL fl99(il) . To convert these ex- 
tinctions to th e Spit zer ban ds we fol l owed the calibrations by 
iFlahertv et all (120071) and ICarpenteri (l200ll) . The relation be- 
tween reddening in 2MASS K s and IRAC bands provided by 
IFlahertv et al. (2007) for p Ophiuchi indicates that no reddening 
is present in the [5.8] - [8.0] color index. Therefore, dereddening 
shifts the objects vertically downwards in Fig. [8] The right hand 
panel of Fig|8] shows the dereddened IRAC colors. Evidently, 
most objects with [3.6] - [4.5] > 0.7 in the left hand side of 
Fig- 18] fall in the "canonical" IRAC Class II area after deredden- 
ing. Therefore, the majority of objects with apparently protostel- 
lar IRAC colors are probably strongly reddened Class II sources. 
We note also that the ISOCAM sub-sample without X-ray coun- 
terparts has the same color distribution as the ISOCAM/DROXO 
sample. This indicates that no bias is present against X-ray prop- 
erties. 

iPrisinzano et al.l (120081) found that Class I protostars in the 
Orion Nebula can be separated into two distinct subclasses: the 
first, la, with a rising SED from K to 8/^m and lower X-ray emis- 



sion level than the second one, lb, characterized by rising SED 
up to 4.5yum. This latter group shows IRAC colors more similar 
to those of Class II objects. The first subclass, Class la, is in- 
stead populated by more embedded protostars with lower and, 
perhaps, more absorbed X-ray emission. They find 23 Class 0- 
Ia YSOs, 22 Class lb YSOs and 148 Class II YSOs. In COUP 
the fraction of Class la on Class II YSOs is ~ 16%, the fraction 
of Class I versus Class II YSOs in Orion is ~ 30%. In DROXO 
the fraction of Class I objects with [3.6]-[4.5] > 1 (after dered- 
dening) is ~ 7% with respect to Class II objects. This sample 
should be com posed by very embedd ed objects similar to Class 
la defined by Prisinzano et al.l (120081) . The number of very em- 
bedded protostars in Core F of p Ophiuchi is thus very low with 
res pect to Class II o b jects . The same evidence has been obtained 
by lj0rgensen et al.l d2008h . In that paper the fraction of Class I 
to Class II is reported to be lower in p Ophiuchi than that present 
in other star-forming regions like Perseus (10% in p Ophiuchi 
vs 90 % in Perseus, respectively), hinting for a different age or 
formation time-scales for these two star-forming regions. 



Three objects (GY92, GY289, GY203) which have been 
classified as Class III by BoOl are found with Spitzer colors sim- 
ilar to Class II YSOs. Furthermore, we find three objects clas- 
sified as Class II YSOs (GY 146, GY 240, GY 450) that are 
bluer than 0.4 mag in ([5.8] - [8.0]) IRAC color . These objects 
may have transition disks (see iKim et al.L 12009, and references 
therein). The star with the reddest IRAC colors located in the top 
right corner of Fig.[8]is WL22/GY174. It is classified as a Class 
II object and proba bly it suffers o f stron g foreground extinction 
as pointed out bv IWilking et al.l (l200ll) . This is supported also 
from our spectral fit to its X-ray spectrum: we find a Nh absorp- 
tion value of 1 .3 • 10 23 cirT 2 , a factor of 5 higher than the average 
of the DROXO sample. 
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6. Results 

6.1. Coronal temperatures and absorption 

Table IA.2I lists the temperatures obtained through a thermal 
model fit to the X-ray spectra as described in Sect. I3.2.2I Most 
of the spectra are reasonably well described by a model with a 
single thermal component; usually we find temperatures higher 
than 1.5 keV and absorption column Nh higher than 10 22 cm -2 . 
In a few cases two or three thermal components are needed to fit 
the spectra, depending on the characteristics of the spectrum and 
its count statistics. For these cases we calculated the average of 
the two or three plasma temperatures weighted by their emission 
measures to obtain a representative mean temperature for these 
spectra. The median of the representative plasma temperatures is 
-3.1 keV with a lcr range between 1.6 and 8 keV (the 10%-90% 
quantile range is 0.8-13 keV). The 10%-90% quantile range of 
the Nh column is 2.8 ■ 10 21 - 7.2 • 10 22 cm -2 and the median is 
2.3 • 10 22 cm" 2 . 

6.2. X-ray luminosities and stellar parameters 

We examine the relation between X-ray luminosity and X-ray to 
bolometric luminosity ratio and mass and effective temperature 
for all X-ray sources for which stellar parameters could be de- 
termined in Sect. 13.31 (Fig. l9l. The X-ray luminosities for the 
sample with known stellar parameters increase from ~ 2 • 10 28 
ergs s- 1 up to ~ 7-10 30 ergs s" 1 in the 3000-5000 K range (Fig.[9] 
top left panel) or ~ 0.1 - 0.7 M (left bottom panel). Most of the 
YSOs have effective temperatures around 3000-5000 K (K type 
stars), with the exception of the three hot stars WL 5, WL 16 and 
WL 19. The fraction of X-ray detections among ISOCAM YSOs 
increases from ~ O.1M to ~ 1M . The boxes in the bottom left 
panel of Fig. [9] are the 10-90% quantile s obtained with ASURV 
software ( Sch mitl [19851 lLavallev et all 1 1 992l) for log L x values 
in four ranges of mass (< 0.2, 0.2 - 0.5, 0.5 - 1.0, > 1.0 M ). 
For these mass ranges the median values of logLx are <28.3, 
29.0, 30.0 and 30.0, respectively. The COUP sample contains 
very few upper limits in each mass range whereas in DROXO 
the fraction of upper limits is ~ 30 % reducing the L x medians 
in each mass range. As discussed below in Sect. 16.31 we sug- 
gest that a fraction of Class III YSOS are likely spurious cloud 
members. We fitted the relation between X-ray luminosity and 
mass excluding these suspect members in the Class III sample 
obtaining logL x ~ 29.9 + 0.2 + l gM( 145±a25) with the same 
procedure used by Preib isch et al.l d2005l) for COUP sample. We 
also fitted X-ray luminosity and mass relation for Class II and 
Class III samples separately, finding similar slopes but different 
normalizations. The relations are 
logLx ~ 29.8 + 0.2 + logM (1 5±03) and 
logLx ~ 30.0 ± 0.3 + logM (L2±a5) 

for Class II and Class III YSOs respectively. The slopes of these 
relations are very similar to t hose f ound for ONC in COUP pro- 
gram (1.44, [Preibisch etal., 2005). The lower normalizations 
that we find in the p Ophiuchi sample with respect to ONC 
(up to 0.6 dex for DROXO Class II objects with respect to 
COUP) can be due to systematic effects like distance estimate 
(for Orion it was used a value of 450 pc, mo re recent estimates 
place ONC to -400 pc, see Sect. 8.1.1 in iMavne & Navloil 
2008, and references therein) and correction for absorption, but 
an intrinsic difference could be present. Thus PMS stars in 
p Ophiuchi are on average less luminous than PMS stars in 
Orion. Because there are upper limits in DROXO, this difference 
should be more marked. Also in the Taurus Molecular Cloud 



(TMO lTelleschT et al.l (120071) find a power law index of — 1.5 
between log Lx and mass. 

The logLx/Lboi ratio increases with decreasing stellar 
masses (Fig. [9j bottom and top right panels) saturating at - 
-3.5 for stars cooler than 5000 K and with masses of - 0.7 
M , slightly below the "canonical" value of -3 observed in 
MS young stars, bu t very similar to the value reported by 
iPreibisch et all (l2005h for Orion YSOs (-3.6). 

From the relation between mass and X-ray luminosity and 
the hypothesis that the r atio logLx/Lboi for PM S stars is satu- 
rated at a level of -3.5, Tellesc hi et alj d2007l) derived an em- 
pirical mass - bolometric luminosity relation for PMS stars 
logLboi - log(M/M ) 149 which is shallower than the relation 
that yields for MS stars. By considering that the saturation limit 
of logLx/Lboi for our sample of YSOs of p Ophiuchi is similar 

to that of TMC and ONC ( 3.5, see fig. [3 and that the slope of 

the relation between mass and L x is quite similar for p Ophiuchi, 
ONC and TMC, it i s suggested that the r elation between mass 
and Lboi found from Telleschi et al] d2007l) for coeval PMS stars 
should apply also for our sample. 

Among the three hottest stars, characterized by effective 
temperatures around 10000 K, we detected the Class III YSOs 
WL 19 and WL 5, while WL 16 (Class II) remains undetected. 
X-ray spectra of WL 19 and WL 5 show quite hot plasma tem- 
peratures (kT= 3.7 keV and 4.5 keV, respectively) and high 
absorption (N# — 10 23 and 6.5T0 22 cm" 2 , respectively). The 
undetected Class II YSO WL 16 is a peculiar object: it con- 
sists of a massive star (L= 25OL , M - 4M ) that illumi- 
nates a circumstellar disk visible only at mid-IR wavelengths 
dRessler & BarsonvLl2003h . It suffers of strong absorption (Ay > 
30) likely due to a foreground screen of cloud material. Taking 
into account this high absorption, the upper limit to luminosity 
is inversely correlated with the plasma temperature: for kT = 
4 keV we obtain log L x £ 28.8, for kT = 1 keV logL x £ 29.6, 
for kT = 0.5 keV logLx ^ 30.7. According to their posi- 
tion in the HR diagram these three objects are intermediate- 
mass pre-MS stars and their X-ray luminosities and L x /Lboi ra- 
tios are in the range typically observed for Herbig Ae/Be stars 
dStelzeretalll2009l) . 

The coole s t obj ect is the binary system WL2/GY128 
dBarsonv et all 120051) . A discrepancy between its spectral type 
and the effective temperature is present in the literature. While 
its spectral type is c o mprise d between K and M as reported 
by [E uhman & Riekd (Q999), its temperature, estimated by 
Natt a et all d2006l) . is very low (-2300 K). As discussed in Sect. 
14.11 WL 2 has undergone a huge flare during DROXO. Its qui- 
escent X-ray luminosity is - 10 29 erg s" 1 which is typical for 
young K-M type stars but unexpectedly high for a low mass 
brown dwarf. Likely the photospheric temperature of this object 
is more similar to that of late K or M-type stars. 

Only two bona fide brown dwarfs are in the field of view. 
We detected GY 310 ( log L x = 28.92 erg s" 1 ), but not GY 141, 
which was detected bv lOzawa et al.l d2005l) during a flare with a 
flux 90 times higher than in a previous Chandra observation. 

6.3. X-ray emission of different YSO classes. 

We examined the X-ray detection rates for YSOs in different 
evolutionary states referring to the YSO classification of BoOl. 
Their catalog comprises 16 Class I, 123 Class II, and 38 bona- 
fide Class III. The latter are classified as Class III YSOs on the 
basis of absence of IR excess and detection in X-ray images 
(ROSAT) and/or radio band (VLA). Given the low ROSAT sen- 
sitivity at Lx - 3 • 10 29 erg s~' the sample of Class III YSOs is 
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Fig. 9. Top left panel: X-ray luminosities (0.3-10. keV) versus effective temperatures of YSOs identified with ISOCAM objects 
(BoOl). Symbols are: squares Class II YSOs, triangles Class III YSOs. Upper limits to luminosities were calculated for ISOCAM 
objects not detected in DROXO (vertical arrows). Top right panel: Lx/Lb i ratios vs T e ff. The upper limit with highest T e ff and 
mass is relative to WL 16 and calculated taking into account an Ay ~ 30 mag and a kT — 1 keV (see text). The detection point 
at the lowest T e s is from WL 2 that exhibited the large flare shown in Fig. [7] Bottom left and right panels: X-ray luminosities and 
Lx/Lboi vs. mass (in units of M ). The boxes in the left panel are the 10%-90% quantile ranges of luminosities derived from the 
Kaplan-Meier estimator of XLFs in four mass ranges (< 0.2, 0.2-0.5, 0.5-1.0, > l.OM ). Horizontal segments mark the median 
of Lx in each mass interval, in the lowest mass bin the median is an upper limit. The solid black line is the relation between mass 
and luminosity fitted between 0. 1 and 2 M , blue and green lines ar e the best fits for Class III and Class II YSOs respectively. The 
dashed gray line is the analogous obtained by Preibisch et al. (2005) to Orion young stars. Horizontal lines in the bottom right panel 
are the values of saturation of log Lx/^boi equal to -3 and -3.5. 



presumably incomplete. Bontemps et al] (1200 ll) report also a list 
of 39 candidate Class III stars which are objects with stellar-like 
colors. They are found in the dense part of the cl oud and have no 
X-ray detection at the ROSAT sensitivity limit. Bontemps et alJ 
d2001l) conclude that these objects are in excess with respect to 
the number of expected field stars because of the high extinc- 
tion in the region where they lie, and for this reason these ob- 
jects were assumed Class III star candidates. For the detection 
rates we considered only the fraction of p Oph members that are 
within the DROXO field-of-view (col. 2 of TableQ]). The number 
of X-ray-detected stars from each YSO group is given in Col. 3, 
and the detection fraction is summarized in Col. 4. The last two 
columns of Table [T]provide the medians of the plasma tempera- 
tures and of X-ray luminosities for each sub-sample. These last 
values were obtained with ASURV and take into account the up- 
per limits for non-detected objects. Before proceeding to a dis- 



cussion of the numbers given in Table[T]an assessment of biases 
and completeness is in order. The ISOCAM sample is not well 
suited to select Class III sources. The evidently low X-ray de- 
tection rate of candidate Class III stars from BoOl suggests that 
most of them are spurious members. We detected six Class III 
candidates with log Lx > 28.3, and found that they have Spitzer 
photometry consistent with Class III objects^. There remain 15 
out of 21 Class III candidates undetected at log Lx < 28.3. Their 
masses are comprised between 0.1 and 0.5 M Q , and their dered- 
dened IRAC colors [3.6] - [4.5] are comprised between 0.0 and 
0.2 mag, whereas the six X-ray-detected Class III candidates 
have [3.6] - [4.5] colors up to 0.4 mag. In the following dis- 



8 Four of six are classifie d as 'YSO' and two of s ix are classifie d as 
'star' i n the Spitzer cat alog, Evans et al. (see 2003b); Padge ttet alJ (see 
2008): lEvanset al] (see llOOgf) for this scheme of classification 
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Table 1. X-ray detection rates and median of logarithm of X- 
ray luminosities and plasma temperatures for the different YSO 
groups following the classification given by Bontemps et al. 
(2001). 



YSO Class 



detection 
fraction 



log L x , med median 
[erg/s] keV 



I 
II 

bona-fide III 
candidate III 



11 

48 
16 
21 



9 
31 
15 
6 



82% 
77% 
94% 
29% 



29.6 
28.8 
29.7 
< 28.3 



4.4 
3.1 
2.4 



Total 



96 



61 



64% 



ferent from a typical coronal singly-absorbed multi-temperature 
plasma. This object has no counterpart in the ISOCAM survey, 
although in the literature is identified with GY 266 and isclas- 
sified as a variable star by Alves de Oliveira & Casali (2008). 
Furthermore, it has an IR counterpart in 2MASS and Spitzer cat- 
alo gs. From IRAC and MIPS photo metry this object is cataloged 
kls a normal star (Evans et al., 2003b). The DROXO spectrum of 
thi s object is shown in Fig. [TT] We modeled the spectrum with 
two thermal components, which were differently absorbed. The 
soft component has kT = 0.47|}g keV (~ 5.5 MK) and is ab- 
sorbed by a Nh column of ~ 3 ( 3 , 4 • 10 20 crrT 2 . This is one of the 
softest thermal components we have determined in the whole 



cussion we will also consider a sub-sample of Class III objects 
without the 15 X-ray-undetected Class III candidates. 

Previous studies of star-forming regions have shown that for 
each given mass range, Class I and I I YSOs show lo wer X- 
ray luminosities than Class III YSOs (iNeuhauser et al 



1995 



Flaccomio et al., 2003; Preibisch et al., 2005; Telleschi et al 



2007; iPrisinzano et a l.. 2008). Figure [10] (top left panel) shows 
the Kaplan-Meier estimators of X-ray luminosity function 
(XLF) of Class I, II, and III objects following the BoOl classi- 
fication. Upper limits are mostly concentrated below the lowest 
detection values for Class II and III YSOs. We plot also the XLF 
of Class III YSOs (dots and dashed line) after excluding 15 up- 
per limits of Class III YSO candidates from BoOl. In this way 
we try to correct the bias introduced by spurious members that 
very likely contaminate the Class III YSO sample as discussed 
above. This corrected Class III XLF is similar to the XLF of 
Class I YSOs. The corrected Class III XLF shows higher lumi- 
nosity levels with respect to Class II XLF. Two sample tests yield 
a probability of > 99% for the two distributions of being differ- 
ent, thus suggesting that in p Ophiuchi Class III YSOs are more 
luminous than Class II YSOs. We observe also that Class I YSOs 
have similar X-ray luminosity levels compared with Class III ob- 
jects. The high X-ray detection rate and high X-ray luminosities 
among Class I objects are surprising. For comparison, in ONC 
Prisinzano et al. (2008) have found that the X-ray luminosities 
of Class I are lower than those of more evolved Class II YSOs. 
Our X-ray bright Class I YSOs could be explained as an effect 
of different mass distributions in different samples but, given that 
we cannot estimate masses of Class I objects with our method, 
we cannot test this hypothesis. 

The distribution of mean plasma temperatures in Fig.fTOldop 
right panel) shows that on average the emitting plasma in Class I 
YSOs is hotter than in Class II and Class III YSOs. The median 
plasma temperatures for Class I, II and III are reported in Table 
[TJ Two-sample tests give a probability of 99% to reject the null 
hypothesis that the distributions of temperature for Class III and 
Class I YSOs are drawn from the same distribution, the probabil- 
ity is ~98% when comparing Class II and Class I temperatures, 
and ~ 90% when comparing Class II and Class III YSOs tem- 
perature distributions. As expected, the distributions of A^h show 
that the absorption is higher on average in Class I YSOs than in 
Class II and III (see Fig. [TOl bottom left and right panels), indi- 
cating circumstellar matter in Class I objects. 

6.4. Source with a soft X-ray excess 

We discuss the case of the source nr. 61, which shows an ex- 
cess of soft photons in the 0.3-1.0 keV band significantly dif- 



sample of DROXO spectra. The hot component has kT = 1.8 2 \ 



2 5 3 - 2 - 



10 22 cm 2 . The emission measure 



keV absorbed by Nh 

ratio of soft to hot component is ~ 1%. The hot component has 
the typical temperature and column density of a YSO in p Oph, 
whereas the weakly absorbed component is unusually soft. 

Idudel et al] (120071) have reported similar examples of X-ray 
spectra with soft excess from the XEST survey. They modeled 
them with a 2-T plasma with individual absorptions and hy- 
pothesized that the soft c omponent be asso ciated to the emis- 
sion from shocks in jets. iBonito et al] (120071) have studied the 
soft X-ray emission that could arise from bipolar jets in YSOs 
through extensive MHD simulations. They find that a jet less 
dense than the ambient medium in which it propagates can emit 
soft X-rays with plasma temperatures of 2-3 MK, very similar 
to that of t he soft component we found. Moreover, Ball v et all 
(2003]) and lFavata et al] (120061) observed that the X-ray emission 
of HH 154 jet is located at the base of the shock and that mor- 
phological changes are detected on a time scale of four years. 
The X-ray luminosity of that jet is 4 • 10 29 erg s _1 . In our case 
the unabsorbed luminosity of the soft component in the spec- 
trum of GY 266 is much lower than that of HH 154 (3.4- 10 27 erg 
s ). Given the angular resolution of EPIC we cannot spatially 
resolve the soft component fro m the hot compon ent. Thanks to 
angular resolution of Chandra (Giidel et al ., 2008) resolved sep- 
arated soft emission from the coronal emission in the star DG 
Tau belonging to the Taurus Molecular Cloud. The coincidence 
of the elongated soft X-ray emission with the direction of the op- 
tical jet represents a clear case for the jet scenario in DG Tau. In 
our case the lack of IR excess is difficult to reconcile with the cir- 
cumstellar disk, which is expected along with the jet. Explaining 
the origin of the X-ray soft component of GY 266 requires fur- 
ther investigation of this star. 



7. Discussion and conclusions 

We described the Deep Rho Ophiuchi XMM-Newton 
Observation (DROXO), aimed at exploring at high sensi- 
tivity the X-ray emission of YSOs in the p Ophiuchi core F 
region. We detected 111 X-rays sources, and for 91 of them we 
obtained a model fit to their X-ray spectra. By using optical and 
IR data we estimated the photospheric parameters of most of the 
sample of X-ray detected YSOs. 

We find 10 unidentified sources and a further 6 sources al- 
ready detected in previous X-ray surveys, but without optical/IR 
counterpart. Three of them (src. num. 45, 48 and 95) show light 
curves with some impulsive variability. The spectra of src. num. 
45 and 48 have a good fit with an absorbed power law, while src. 
num. 95 has a spectrum compatible with a thermal model. The 
lack of 2MASS counterparts could be a hint that they are of ex- 
tragalactic nature, but, given the sensitivity of 2MASS catalog, 
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Fig. 10. Top left panel: Kaplan-Meier estimators of cumulative distributions of Lx for Class I YSOs (solid line), Class II YSOs 
(dashed line), Class III YSOs (dotted line) and corrected XLF (dotted line + points) of Class III YSOs from BoOl . Top right panel: 
cumulative distribution of mean plasma temperatures obtained from fit to spectra for Class I (solid line), Class II (dashed line), 
Class III (dotted line) objects. Bottom left and right panels: cumulative distributions of A^h column absorption among the different 
classes (left) and scatter plot of mean plasma temperatures vs. A^h (right). Circles are Class I YSOs, squares are Class II YSOs, and 
triangles are Class III YSOs. 



it is not possible to rule out they are very low mass PMS stars or 
even brown dwarfs, especially src. num. 95. 

The sensitivity of the survey is fx ~ 5 ■ 10~ 15 erg s _1 cm -2 
in flux and log Lx (erg s _1 ) ~ 27.9 in luminosity, but it strongly 
depends on the local absorption, which is largely variable in the 
field-of-view. The sample of 96 classified YSOs given by BoOl 
in the DROXO field of view has allowed us to explore the X-ray 
emission from Class I, II, and III YSOs of p Ophiuchi. The X- 
ray properties of p Ophiuchi PMS stars obtained from DROXO 
were compared with those of the Orion Nebula Cloud PMS stars 
studied in the COUP survey. When fitting the relation between 
mass and X-ray luminosity with a power law, the index that we 
find for DROXO sa mple is quite similar to that found for ONC 
and TMC dPreibisch et aUl2005tlTelleschi et aUl2007l) . With re- 
spect to COUP we find a lower normalization, suggesting that 
on average the PMS stars in the range of masses between 0.1-2 
M G in our sample are less luminous than their analogs in COUP. 
Young stellar objects in p Ophiuchi exhibit a saturation of the 
ratio Lx/Lboi near the "canonical" value of ~ 1CT 3 for masses 
between 0.5 and 1 M Q , while stars with masses below 0.5 Mo 
show a lower limit of saturation Lx/Lboi ~ 10~ 3 5 . This is con- 



sistent with what was found by [Preibisch et alj d2005l) for the 
COUP sample. 

We detected a large fraction of YSOs in the field of view, 
94% of which are bona fide Class III stars, 77% of Class II and 
82% of Class I YSOs, respectively. We confirm the high detec - 
tion rate among Class I YSOs found by Imanishi et al] (1200 ll) . 
The detection rate in our Class I sampl e is higher than t h e ana- 
log rate found in the COUP survey by iPrisinzano et alj {2008) 
(82% vs. 62%) d espite the higher sensiti vity of COUP with re- 
spect to DROXO. IPrisinzano et alJd2008h make a distinction be- 
tween Class la (characterized by a rising SED from K to 8yum) 
and Class lb (rising SED up to 4.5/mi). The detection rates in 
these two subclasses are 44% and 82% respectively. X-ray lu- 
minosities of Class lb objects are higher than those of Class la 
objects. The rate of X-ray detection of Class lb objects in Orion 
is identical to that of our sample of Class I objects in p Ophiuchi. 
With support of Spitzer photometry we suggest that only a small 
fraction of our Class I stars are deeply embedded objects. From 
X-ray data the high detection rate suggests that our Class I sam- 
ple could be mainly formed by Class lb YSOs as defined by 
Prisinzano et al. (2008). These objects in Orion have X-ray lu- 
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Fig. 11. X-ray spectrum of source # 61 (GY 266) with a soft 
component in the 0.3-1 keV energy range, less absorbed than 
the hot component. 



minosities more similar to Class II objects. In p Oph we find that 
our nine detected Class I stars are on average more luminous 
tha n Class II stars. 

Bontemp s et al.1 (1200 ll) list 39 stars defined as Class III can- 
didates, 21 of which are surveyed in DROXO. Six out of 21 are 
detected in DROXO, providing support for their PMS nature. 
The other 15 are undetected below logLx = 28.3. Likely they 
are not PMS stars and thus not members of the p Oph cloud. By 
means of XLFs we evaluated the level of emission of Class I, II, 
and III YSOs. After excluding the suspect Class III YSOs, we 
find that Class I YSOs emit at the same level of Class III YSOs 
(median log Lx ~ 29.6 and 29.7, respectively) and both samples 
are more luminous than Class II YSOs (median log Lx < 28.8). 
In COUP Class III objects are more luminous than Class I and II 
for masses between 0.5 and 1.2 M0, while for masses in 0.1-0.5 
MO Class lb source s are slightly more lum inous than Class II 
and III (see Fig. 8 in lPrisinzano et aUl2008h . 

The analysis of X-ray spectra leads us to conclude that the 
mean plasma temperatures of Class I protostars are higher than 
in Class II and Class III YSOs. We find a clear trend of de- 
creasing plasma temperatures passing from Class I to Class 
III objects, differently from what is observed in ONC where 
iPrisinzano et alJ ([2008) do not find a significant evolution of 
temperatures from Class I to Class III objects. Four Class III 
stars and one Class II star (GY 1 12, SR 12A, GY 296, GY 380 
and GY 3) have spectra noticeably softer than the rest of YSOs 
with mean temperatures around 0.7-0.8 keV. 

The absorption derived from X-ray spectra is higher in 
Class I than in Class II and III. Given the high absorption among 
Class I YSOs, it is impossible to determine if a soft thermal com- 
ponent in their spectra is completely absorbed or is not present at 
all. This soft component would reduce the intrinsic kT for Class I 
objects, making them more similar to the ones in the ONC. 
However, higher absorption should lead to an underestimate of 
the luminosities of heavily absorbed spectra like those of Class I 
objects, while we observe a luminosity for them higher than that 
of in Class II and Class III YSOs. Furthermore, Prisinzano et al. 
(2008) have found in extensive simulations that varying Nh does 
not significantly influence the XLF 



The star GY 266 shows a peculiar X-ray spectrum composed 
by two thermal components that are differently absorbed. The 
hot and heavily absorbed component is similar to that found in 
other PMS stars, while the soft one, less absorbed, could arise 
from unresolved jet shocks. 
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Appendix A: Online tables 



Table A.l. List of sources detected in the DROXO EPIC images. The columns refer to X-ray positions, position accuracy, summed exposure times 
from the three EPIC c ameras, count rates scal ed to MOS detector units (see Sect. |3. It with l<x errors, 2MASS identifier, ISOCAM identifier, SED 
classification given by Bontemps et al. ( 2001), literature names, flag column for identifications (see footnote). Coordinate errors take into account 
position uncertainties calculated by the detection code and the best-match radius from the optical/IR objects. Class III candidates from BoOl are 
indicated by "III?". 
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I 


36 


16:27:09.4 


-24:43:20.3 


1.6 


3.6 


259.4 


3.8 ±0.1 


16270931-2443196 








S 


37 


16:27:09.4 


-24:40:22.6 


2.8 


1.7 


592.4 


0.37 ± 0.04 


16270935-2440224 


107 


II 


GY213 


I 


38 


16:27:09.5 


-24:37:19.8 


1.2 


3.2 


455.1 


10.9 ± 0.2 


16270943-2437187 


108 


I 


EL29/GY214 


I 


39 


16:27:11.2 


-24:40:47.0 


1.7 


1.4 


446.0 


2.1 ±0.09 


16271117-2440466 


112 


II 


GY224 


I 


40 


16:27:11.7 


-24:38:32.9 


1.5 


2.0 


572.4 


3.18 ±0.09 


16271171-2438320 


114 


III 


WL19/GY227 


I 


41 


16:27:12.7 


-24:40:52.9 


2.9 


1.2 


611.5 


0.084 ± 0.02 










X 


42 


16:27:13.9 


-24:43:33.6 


3.0 


3.5 


592.5 


0.29 ± 0.04 


16271382-2443316 


117 


II 


GY235 


I 


43 


16:27:15.1 


-24:51:39.3 


1.2 


11.5 


324.5 


54.4 ± 0.5 


16271513-2451388 






WSB 46 


S 


44 


16:27:15.5 


-24:30:53.6 


1.8 


9.2 


269.8 


0.88 ± 0.08 


16271551-2430536 


119 


II 


IRS35/GY238 


I 



Table A.l. continued. 



DROXO source number 


R.A. 


DEC 


Pos. Error 


Off-axis 


Exp. Time 


Rate 


2MASS id. 


ISOCAM id.° 


CLASS 


Name 


Flag 6 




J2000 


J2000 


n 




ks 


ct/ks 










45 


16:27:15.7 


-24:33:04.0 


4.3 


7.1 


329.9 


0.64 ± 0.08 




- 






X 


46 


16:27:15.9 


-24:38:43.7 


1.1 


1.4 


597.1 


13.5+0.2 


16271569-2438434 


121 


II 


WL20/GY240 


I 


47 


16:27:16.4 


-24:31:15.8 


2.4 


8.9 


281.2 


3.8 ± 0.2 


16271643-2431145 


- 






S 


48 


16:27:17.4 


-24:36:25.6 


4.7 


3.7 


523.3 


0.35 ± 0.05 




- 






X 


49 


16:27:18.0 


-24:28:53.4 


2.0 


11.2 


63.7 


60.9 ± 2. 


16271817-2428526 


125 


III 


WL5/GY246 


I 


50 


16:27:18.4 


-24:54:55.2 


2.3 


14.8 


225.9 


6.1 ±0.2 


16271836-2454537 


- 






S 


51 


16:27:18.5 


-24:39:15.4 


1.8 


1.0 


616.5 


1.4 ± 0.06 


16271838-2439146 


127 


I 


GY245 


I 


52 


16:27:18.5 


-24:29:07.7 


2.3 


11.0 


58.4 


1.6 ±0.4 


16271848-2429059 


128-129 


II— II 


WL3/GY249-WL4/GY247 


I 


53 


16:27:19.6 


-24:41:41.4 


0.8 


1.7 


628.2 


100.0 ± 0.5 


16271951-2441403 


130 


III 


SR12/GY250 


I 


54 


16:27:21.5 


-24:41:43.6 


1.3 


1.9 


609.0 


0.6 ± 0.05 


16272146-2441430 


132 


II 


IRS42/GY252 


I 


55 


16:27:21.9 


-24:43:37.0 


1.1 


3.7 


615.8 


11.3 ±0.2 


16272183-2443356 


133 


III 


GY253 


I HH 


56 


16:27:22.1 


-24:29:54.1 


1.9 


10.3 


232.9 


2.8 ±0.1 


16272180-2429533 


134 


I 


WL6/GY254 


I V 


57 


16:27:23.0 


-24:48:08.5 


2.6 


8.1 


454.9 


1.8 ± 0.09 


16272297-2448071 


- 






ss 


58 


16:27:24.4 


-24:50:56.4 


6.6 


10.9 


342.7 


0.5 ± 0.06 




- 






ug 


59 


16:27:24.7 


-24:29:39.0 


2.5 


10.6 


219.5 


5.0 ± 0.2 


16272463-2429353 


- 






s o. 

o n 


60 


16:27:26.5 


-24:39:23.3 


1.0 


2.3 


614.9 


4.08 ±0.1 


16272648-2439230 


140 


11 


GY262 




61 


16:27:27.0 


-24:32:17.4 


2.7 


8.2 


274.5 


2.7 ±0.1 


16272706-2432175 


- 




GY266 


s : - 


62 


16:27:27.1 


-24:40:51.9 


0.9 


2.5 


644.4 


22.0 ± 0.2 


16272693-2440508 


141 


I 


IRS43/GY265 


T * 

1 CD 


63 


16:27:27.4 


-24:31:17.1 


2.1 


9.2 


271.3 


5.6 ± 0.2 


16272738-2431165 


142 


II 


VSSG25/GY267 


c« 

I & 


64 


16:27:28.1 


-24:39:34.7 


1.1 


2.6 


617.5 


13.7 ±0.2 


16272802-2439335 


143 


I 


IRS44/GY269 


:s r 


65 


16:27:28.3 


-24:27:21.3 


3.2 


13.0 


57.5 


1.5 ±0.2 


16272844-2427210 


144 


II 


IRS45/GY273 


si 


66 


16:27:28.5 


-24:54:34.3 


2.5 


14.7 


236.2 


3.9 ± 0.2 


16272873-2454317 


- 






67 


16:27:29.5 


-24:39:17.0 


1.4 


3.0 


612.4 


0.36 ± 0.03 


16272943-2439161 


145 


I 


IRS46/GY274 


'§ 


68 


16:27:30.0 


-24:27:42.8 


2.8 


12.8 


20.2 


0.9 ± 0.4 


16272960-2427419 


- 






I o 


69 


16:27:30.0 


-24:33:37.2 


1.7 


7.2 


353.6 


2.9 ±0.1 


16272996-2433365 


146 


III? 


GY278 


70 


16:27:30.6 


-24:52:09.7 


3.7 


12.4 


315.0 


3.3 ±0.1 




- 






S 3 


71 


16:27:30.9 


-24:47:27.7 


1.1 


8.0 


446.8 


12.7 ±0.2 


16273084-2447268 


149 


Ill 




72 


16:27:31.2 


-24:34:04.3 


2.0 


6.9 


366.2 


1.7 ± 0.09 


16273105-2434032 


148 


III? 


GY283 


lin 


73 


16:27:31.3 


-24:31:11.1 


2.3 


9.5 


283.2 


0.52 ± 0.06 




- 






si 


74 


16:27:32.6 


-24:44:59.6 


3.6 


6.0 


561.5 


0.21 ± 0.04 


16273272-2445004 


- 






75 


16:27:32.8 


-24:33:24.8 


1.7 


7.6 


343.0 


3.2 ±0.1 


16273267-2433239 


152 


Ill 


GY289 


I f? 


76 


16:27:32.9 


-24:32:35.5 


2.2 


8.4 


318.9 


4.7 ± 0.2 


16273285-2432348 


154 


II 


GY291 


is- 


77 


16:27:33.4 


-24:41:15.8 


1.0 


3.9 


477.4 


26.4 ± 0.3 


16273311-2441152 


155 


II 


GY292 




78 


16:27:35.6 


-24:38:34.1 


2.0 


4.5 


401.8 


1.5 ± 0.09 


16273526-2438334 


156 


III? 


GY295 


I Uj 


79 


16:27:35.8 


-24:45:34.1 


1.8 


6.9 


519.2 


1.5 ± 0.07 


16273566-2445325 


157 


Ill 


GY296 


I 


80 


16:27:37.4 


-24:42:39.6 


1.9 


5.3 


443.8 


1.5 ±0.08 


16273724-2442380 


161 


II 


GY301 


I 


81 


16:27:38.2 


-24:30:43.1 


3.2 


10.6 


263.6 


1.7 ±0.1 


16273812-2430429 


- 




IRS50/GY306 


S 


82 


16:27:38.3 


-24:37:01.6 


1.1 


5.8 


514.1 


11.8 ±0.2 


16273832-2436585 


163 


11 


IRS49/GY308 


I 


83 


16:27:38.7 


-24:38:39.5 


1.7 


5.2 


543.5 


0.64 ± 0.05 


16273863-2438391 


164 


II 


GY310 


I 


84 


16:27:39.0 


-24:40:19.6 


1.8 


5.0 


302.6 


0.73 ± 0.07 


16273894-2440206 


165 


II 


GY312 


I 


85 


16:27:39.1 


-24:47:21.6 


3.2 


8.8 


403.6 


3.6 ±0.1 




- 






S 


86 


16:27:39.3 


-24:39:14.9 


0.9 


5.2 


396.4 


43.0 ± 0.4 


16273942-2439155 


166 


II 


GY314 


I 


87 


16:27:39.9 


-24:43:13.6 


1.3 


6.1 


426.8 


8.84 ± 0.2 


16273982-2443150 


167 


I 


IRS51/GY315 


I 


88 


16:27:40.6 


-24:46:55.3 


3.5 


8.7 


444.5 


0.23 ± 0.03 










u 


89 


16:27:41.6 


-24:35:37.8 


4.0 


7.2 


362.8 


0.64 ± 0.07 


16274149-2435376 


169 


III? 


GY322 


I 


90 


16:27:42.8 


-24:38:51.6 


2.5 


6.0 


521.6 


0.63 ± 0.05 


16274270-2438506 


172 


II 


GY326 


I 


91 


16:27:43.7 


-24:31:27.7 


2.9 


10.6 


269.3 


2.7 ±0.1 










S 


92 


16:27:44.4 


-24:48:57.4 


5.1 


10.8 


319.9 


1.1 ±0.08 










S 



Table A.l. continued. 



DROXO source number 


R.A. 
J2000 


DEC 
J2000 


Pos. Error 


Off-axis 
' 


Exp. Time 
ks 


Rate 
ct/ks 


2MASS id. 


ISOCAM id." 


CLASS 


Name 


Flag 6 


93 


16:27:47.1 


-24:45:35.4 


1.5 


8.8 


430.2 


3.77+0.1 


16274709-2445350 


177 


II 


GY352 


I 


94 


16:27:49.2 


-24:39:07.2 


3.2 


7.4 


448.3 


0.81 ±0.06 




- 






S 


95 


16:27:50.3 


-24:31:47.2 


4.6 


11.3 


258.8 


1.1+0.1 




- 






X 


96 


16:27:51.7 


-24:47:45.3 


7.0 


11.0 


382.7 


0.41 ± 0.06 




- 






u 


97 


16:27:51.9 


-24:31:45.4 


2.4 


11.6 


252.8 


0.41 + 0.06 


16275180-2431455 


182 


I 


IRS54/GY378 


I 


98 


16:27:52.0 


-24:46:30.2 


1.4 


10.2 


409.9 


6.18 ±0.1 


16275191-2446296 


183 


III 


GY377 


I 


99 


16:27:52.2 


-24:40:51.3 


1.1 


8.1 


284.0 


158 ± 1 


16275209-2440503 


184 


III 


IRS55/GY380 


I 


100 


16:27:55.8 


-24:44:51.1 


3.5 


10.1 


417.9 


0.88 ± 0.08 


16275565-2444509 


186 


III? 


GY398 


[ 


101 


16:27:57.9 


-24:40:02.4 


1.4 


9.3 


418.9 


5.99 ±0.1 


16275782-2440017 


188 


Ill 


GY410 




102 


16:27:60.0 


-24:48:20.8 


1.3 


12.8 


327.6 


8.51 ±0.2 


16275996-2448193 


- 






s 


103 


16:27:60.0 


-24:42:41.5 


8.3 


10.1 


375.2 


0.8 ± 0.08 




- 








104 


16:28:02.0 


-24:49:55.7 


6.5 


14.2 


263.8 


0.7 ± 0.09 




- 






S hrt 


105 


16:28:04.7 


-24:34:55.8 


1.6 


12.1 


306.7 


16.9 ±0.3 


16280464-2434560 


191 


III? 


GY463 


1 1 


106 


16:28:04.8 


-24:37:10.5 


2.3 


11.3 


301.1 


3.9 ±0.1 


16280478-2437100 








S 2 
<,3. 


107 


16:28:05.0 


-24:50:18.7 


9.1 


14.9 


51.4 


1.0 ± 0.2 










108 


16:28:07.0 


-24:48:29.9 


3.9 


14.1 


127.1 


1.4 ±0.1 










s 

s : - 


109 


16:28:08.0 


-24:45:10.9 


3.5 


12.7 


235.9 


0.48 ± 0.08 


16280810-2445121 








110 


16:28:08.8 


-24:48:41.5 


5.5 


14.6 


122.9 


0.58 ±0.1 










111 


16:28:10.0 


-24:39:31.7 


5.5 


12.1 


328.9 


0.28 ± 0.05 










SE. 



Notes. 

" Identifier from dBontemps et aUl200lh . 

.1 indicates a match in the ISOCAM survey BoOl. IRS 50 (Src. num 81) is liste d i n BoOl, but undetec ted in the ISOCAM survey. S indicates a match with one or more Spitzer objects. X indicates 
a match only with X-ray sources detected in the surveys Flaccomio et al. (2006) or Ozawa et al. U indicates a source without any X-ray, optical and IR counterpart. 
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Table A.2. X-ray Parameters from model fit to source spectra. The errors are quoted at the 90% confidence region of best fit. Unabsorbed X-ray 
fluxes and luminosities are calculated in the 0.3-8 keV energy band. 



DROXO Src. num. CCDs 


Instr 


SNR Model N H 


kT 1 


log E. M. 1 kT 2 


log E. M. 2 


kT3 


E. M. 3 


log F x log L x x 1 


D.o.F. Prob 






10 22 cm" 2 


keV 


cirT 3 keV 


cm 4 


keV 


cirT 3 


erg s -1 cirT 2 erg s -1 





77" 

'0.63 

7 13 

"•'03 

3 3 3 5 

3.1 

4 5 50 

> 10 

Q 146.0 
y - 1 4.9 

> 10 

0.7 6 ' 4 
5¥- 76 

14 

> 10 
1.6" 
0.5'J 

2 2 63 

^0.9 

3 041 

3 l 74 

2 6 28 



5L67 |f 

50.SS 5 ' 00 
52 - 23 iis 

51 - 78 1 

52. 2 1 52 - 27 
51 66^ 

5169 

51.36P 4 
51 45 5155 

51.36 

51.94g 326 
5 3 3 9 5345 

JJ,J7 S! .33 

51 37 5L86 

J1-J '51.08 

50-88^ 
51.84 2 : >> 

51.11512 

52 - 94 l'l 

52 79" 2 ^ 

52.63|-| 
52 52" " 
51.48^; f 
51 61 5192 



1 


M2PN 


ALL 


L 


IT 


2 


M1PN 


ALL 


L 


IT 


3 


M1M2PN 


PN 


H 


IT 


4 


M1M2PN 


PN 


H 


IT 


5 


M1M2PN 


ALL 


L 


IT 


6 


M1M2PN 


ALL 


L 


IT 


8 


M1M2PN 


ALL 


L 


IT 


9 


M1M2PN 


ALL 


H 


2T 


10 


M1M2PN 


ALL 


H 


IT 


11 


M2PN 


ALL 


H 


IT 


12 


M2PN 


ALL 


L 


IT 


13 


M1M2PN 


ALL 


H 


IT 


14 


M1M2PN 


ALL 


H 


IT 


15 


PN 


PN 


H 


IT 


17 


M1M2PN 


ALL 


L 


IT 


19 


M2PN 


PN 


L 


IT 


22 


M1M2PN 


ALL 


L 


IT 


23 


M2PN 


PN 


L 


IT 


25 


M1M2PN 


ALL 


H 


IT 


27 


M1M2PN 


M1PN 


H 


IT 


28 


M1M2PN 


ALL 


L 


IT 


30 


M1M2PN 


PN 


H 


IT 


31 


M1M2PN 


ALL 


H 


IT 


33 


M1M2 


ALL 


H 


IT 


34 


M1M2PN 


ALL 


H 


IT 


35 


M2PN 


PN 


H 


IT 


36 


M1M2 


ALL 


L 


IT 


37 


M1M2PN 


ALL 


L 


IT 


38 


M1M2PN 


M2 


H 


IT 


39 


M1M2PN 


M1M2 


H 


IT 


40 


M1M2PN 


ALL 


L 


IT 


43 


M2PN 


PN 


H 


3T 


44 


M2PN 


ALL 


H 


IT 


45 


M2PN 


ALL 


L 


IT 


46 


M1M2PN 


PN 


H 


IT 


47 


M2PN 


ALL 


H 


IT 


49 


M2 


ALL 


H 


IT 


50 


M2PN 


ALL 


H 


IT 


51 


M1M2PN 


ALL 


L 


IT 


52 


M2 


ALL 


H 


IT 


53 


M1M2PN 


PN 


H 


3T 


54 


M1M2PN 


M1PN 


H 


1T+3T 


55 


M1M2PN 


PN 


L 


IT 


56 


M2PN 


ALL 


H 


IT 


57 


M2PN 


ALL 


H 


IT 



24 



1.1' 



o-Ki 

0.93; 97 

o.t 7 „ 
09° 59 
3 9 4 - 8 
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y - L AA 

°- 38 035 

0.44 jj 
7'' 
25' 7 
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.7 
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0-05™ 
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0.4 
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53 



6.6 



15.0 



.66 



3.7 4 ;; 
2.40L I7 
2.8=| 
5.4 
3.01| 

3 9 5 5 

4 5 5 5 

3.8|* 
2 2 6 

■^•"1.6 
410.42 
U - Hi 0.40 

2.5 33 



0.29! 



0.22 



2. 
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1.6* 



^2.0 
33.0 



Mr 



L3_ 



52.4 

52.03 

52.04^ 

52.96? 2jr 

53.2 
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52 24 52 55 
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5 1.8 If " 
53.02" 05 
52.7i a91 
53.65 
51.57 
52.31 
52-96„ 
52.98t° 3 
52.26 5241 
52.98: 
52.81 5272 
51 42 5L5 



■M 

m 

52.09 
,53.15 



— 


— 


— 


— 


1.8e-14 


28.5 


48 


37 


0.11 


— 


_ 


_ 


_ 


3.4e-14 


28.8 


14.5 


10 


0.15 


— 


_ 


_ 


_ 


1.4e-12 


30.4 


262.3 


221 


0.03 


_ 


_ 


_ 


_ 


3.5e-12 


30.8 


132.4 


125 


0.31 


_ 


_ 


_ 


_ 


7.3e-15 


28.1 


1.9 


3 


0.6 


_ 




_ 


_ 


1.8e-13 


29.5 


67.7 


90 


0.96 


_ 


_ 


_ 


_ 


5.8e-14 


29 


9.7 


8 


0.28 


°- 70 o:67 


52- 10g^ 


_ 
_ 


_ 
— 


2e-13 
3.2e-14 


29.5 
28.7 


168.2 
21.3 


144 
12 


0.082 
0.046 


_ 


_ 


— 


_ 


6.4e-14 


29 


11.4 


9 


" 0.25 


_ 


_ 


_ 


_ 


2.3e-14 


28.6 


1.7 


2 


3 0.43 
n 0.065 


_ 


_ 


_ 


_ 


1.9e-14 


28.5 


21.4 


13 


_ 


_ 


_ 


_ 


6e-14 


29 


10.6 


9 


2. 0.3 


_ 


_ 


_ 


_ 


2.5e-12 


30.6 


85.4 


74 


£0.17 

: 0.86 


_ 


_ 


_ 


_ 


1.7e-14 


28.5 


6.2 


11 


_ 


_ 


_ 


_ 


5.1e-15 


27.9 


7.5 


6 


£0.28 


_ 


_ 


_ 


_ 


5.5e-14 


29 


49.1 


34 


! 0.045 


- 


- 


- 


- 


9.8e-15 


28.2 


18.8 


17 


« 0.34 


_ 


_ 


_ 


_ 


6.5e-13 


30 


402.5 


367 


cf 0.098 


_ 


_ 


_ 


_ 


4.8e-13 


29.9 


18.4 


13 


rt - 14 


_ 


_ 


_ 


_ 


2.9e-14 


28.7 


27.4 


20 


§0.12 





_ 







2.9e-13 


29.7 


157.8 


126 


g 0.029 


- 


- 


- 


- 


2.3e-13 


29.6 


182.9 


154 


O 0.056 










2.2e-14 


28.6 


9 


9 


" _ 0.43 


_ 


_ 


_ 


_ 


3.9e-14 


28.8 


8 


9 


S. - 53 


_ 


_ 


_ 


_ 


2.5e-13 


29.6 


34.7 


27 


I? 0.15 
g.0.56 
$f 0.27 


_ 


_ 


_ 


_ 


9.3e-14 


29.2 


91.2 


94 


_ 


_ 


_ 


_ 


7.2e-14 


29.1 


10 


8 


_ 


_ 


_ 


_ 


l.le-12 


30.3 


71.3 


48 


§ 0.016 


_ 


_ 


_ 


_ 


1.6e-13 


29.4 


14.8 


14 


SL0.39 


_ 


_ 


_ 


_ 


4.6e-13 


29.9 


225.8 


209 


^ 0.2 


0-1883 


52 2 1 524 ' 


61° 64 

u ' ul 0.58 


52-57£63 


5.9e-13 
1.3e-13 


30.0 
29.4 


284.3 
18 


245 
9 


^0.02 
0.035 


_ 


_ 


_ 


_ 


6.3e-14 


29 


19.5 


12 


0.078 


_ 


_ 


_ 


_ 


8.3e-13 


30.2 


251.1 


203 


0.012 


- 


- 


- 


- 


4.4e-13 


29.9 


31 


24 


0.15 










4.1e-12 


30.8 


61.4 


59 


0.39 










2.6e-14 


28.7 


21.1 


19 


0.33 










1.8e-13 


29.5 


71.9 


65 


0.26 










6.3e-13 


30 


15 


8 


0.059 


3 5 4 - 3 
3.1 


52 50 52 ' 55 


910.94 
u ' 71 0.88 


52-94«f 


8e-13 
1.3e-13 


30.1 
29.4 


587.8 
29.9 


433 
33 


le-06 
0.62 










7.4e-13 


30.1 


509.2 


531 


0.74 










6.4e-13 


30 


53.2 


37 


0.041 










1.7e-14 


28.5 


6.8 


9 


0.66 



Table A.2. continued. 



DROXO src. num. CCDs 



Instr 



SNR Model 



N H 

10 22 cm- 2 



kT 1 

keV 



log E. M. 1 

cm" 3 



59 


M2PN 


ALL 


H 


IT 


60 


M1M2PN 


M2PN 


H 


IT 


61 


M1M2PN 


PN 


L 


2NH 


62 


M1M2PN 


PN 


L 


IT 


63 


M2PN 


ALL 


H 


IT 


64 


M1M2PN 


PN 


H 


IT 


65 


M2 


ALL 


L 


IT 


66 


M2PN 


M2 


H 


IT 


67 


M1M2PN 


PN 


L 


IT 


68 


M2 


M2 


H 


IT 


69 


M2PN 


ALL 


H 


IT 


71 


M1M2PN 


ALL 


H 


IT 


72 


M2PN 


ALL 


H 


IT 


75 


M2PN 


ALL 


H 


IT 


76 


M2PN 


ALL 


H 


IT 


77 


M1M2PN 


M1M2 


H 


2T 


78 


M1M2PN 


M1M2 


H 


IT 


79 


M1M2PN 


ALL 


L 


IT 


80 


M1M2PN 


M2 


L 


IT 


81 


M2PN 


ALL 


H 


IT 


82 


M1M2PN 


ALL 


H 


IT 


83 


M1M2PN 


M2PN 


H 


IT 


84 


M1M2PN 


ALL 


H 


IT 


85 


M1M2PN 


PN 


H 


IT 


86 


M1M2PN 


M1M2 


H 


2T 


87 


M1M2PN 


Ml 


H 


IT 


88 


M1M2PN 


ALL 


L 


IT 


91 


M2PN 


PN 


L 


IT 


92 


M1M2PN 


ALL 


L 


IT 


93 


M1M2PN 


ALL 


H 


IT 


94 


M1M2PN 


ALL 


L 


IT 


95 


M2PN 


ALL 


L 


IT 


96 


M1M2PN 


ALL 


L 


IT 


97 


M2PN 


ALL 


L 


IT 


98 


M1M2PN 


ALL 


H 


IT 


99 


M1M2PN 


Ml 


H 


2T 


100 


M1M2PN 


ALL 


H 


IT 


101 


M1M2PN 


ALL 


L 


IT 


102 


M1M2PN 


ALL 


H 


2T 


104 


M1M2PN 


ALL 


L 


IT 


105 


M1M2PN 


ALL 


L 


IT 


106 


M1M2PN 


ALL 


L 


IT 


107 


Ml 


ALL 


L 


IT 


108 


M1M2 


ALL 


H 


IT 


109 


M1M2PN 


ALL 


L 


IT 


110 


M1M2 


ALL 


L 


IT 



~Y6 n — 

z..u 2 j 

3 3 34 
2 5 3 - 2 

2 I 

4 4 4 -* 

^4.3 

75 089 

u - '^0.65 

7.2™ 

1 o'- 9 

~ 
6 1" 

1 9 2 - 8 
6 3 74 

5.6 

o.8i 55 8 7 

2 l 2 - 4 

1 9 20 

1-8 
32 045 

U,J Q.23 

1.0' 



5 2 9 - 1 
4.3 4 ;J 
1.8- 
3 2 3 3 
2.71 
5.8g 
4.8 
6° 74 

Urn 

u.u 37 

2 2 41 
5 6 5 

2J 25 

2.2 ' 
3.2 

3 n 

74 86 

2.6 3 " 
07 §:§ 

u -'0.6 

-2 764 

J -'2.5 

2.1 30 

2 7 2:5 
^•'2.6 

3 2 4 - 7 

27 

> 10 

"• 7J 0.68 
3 o4.6 
i -°3.2 

> 10 

> 10 
17 

3 9 4 8 

J ' y 3.2 

> 10 

> 10 
4.2 
2.5 

2 2 24 
0.8tS 



CO -2052.49 
Jt,J7 52.2l 

53 01 53 04 
5 2 3 3 52 5 

53.29|| 

52 17' 224 
53.4=g 
51.61 
51.3 gf 
51-82?- 
5 2 3 52 - 57 
52.64*™ 
52.67|g 
51.# 

5235 il 

53 22^'^ 

JJ HI? 

51.65; 
52.03*" 
52 14 5229 
52.14| 2§ 
52.86 52 ; 89 
51.78^[; 87 
51.88 
52.19g 24 
52.75|| 
52.86 2 ? 4 
51.38f 28 
52.25= 
51.60= 
52.27= 

51.33f 47 
52.49f 88 

53.02|| 
51.14 1:24 



kT 2 


log E. M. 2 


■ [ 3 


log E. M. 3 


log F x 


logL x 


n ' 

X 


D.o.F. 


Prob 


keV 


cm - 


keV 


cm - 


erg s cm 


erg s 








— 


— 


— 


— 


2.4e-13 


29.6 


25.7 


22 


0.26 


— 


— 


— 


— 


9.2e-13 


30.2 


164.6 


165 


0.5 


— 


— 


— 


— 


1.4e-13 


29.4 


41.5 


51 


0.83 


— 


— 


— 


— 


1.6e-12 


30.4 


835.5 


854 


0.67 


— 


— 


— 


— 


l.le-13 


29.3 


48 


34 


0.057 


— 


— 


— 


— 


2.5e-12 


30.6 


393.2 


354 


0.074 


— 


— 


— 


— 


3.8e-14 


28.8 


4.3 


4 


0.36 


— 


— 


— 


— 


1.4e-14 


28.4 


5 


4 


0.29 


— 


— 


— 


— 


6.8e-14 


29.1 


38.3 


25 


0.043 










1.4e-13 


29.4 


0.1 


3 


0.99 


- 


- 


- 


- 


4.2e-13 


29.9 


53.7 


36 


i- 1 0.029 


— 


— 


— 


— 


3.5e-13 


29.8 


220.2 


197 


gO.12 


— 


— 


— 


— 


5.6e-14 


29 


10.2 


14 


g0.75 


— 


— 


— 


— 


1.8e-13 


29.5 


35.2 


30 


CD ^ r. a 

3. 0.24 






— 


— 


1.9e-13 


29.5 


28 


29 


" 0.52 

rr 0.21 

£0.27 

CD 

£0.21 
&■ 0.72 


3.0 3 ; 3 
— 


53.O853 3 

— 


— 
— 


— 
— 


2.1e-12 
3.3e-14 


30.6 
28.8 


270.6 
12.3 


253 
10 


— 


— 


— 


— 


7.7e-14 


29.1 


91.9 


82 










1.2e-13 


29.3 


24.2 


29 


- 


- 


- 


- 


9.5e-14 


29.2 


15.5 


12 


3* 0.21 


— 


— 


— 


— 


5.5e-13 


30 


291.4 


239 


3 0.012 


— 


— 


— 


— 


4.8e-14 


28.9 


17.9 


15 


j»0.27 


— 


— 


— 


— 


8.4e-14 


29.2 


31 


14 


9 0.0055 

C ^ rr 

O 0.55 






— 


— 


1.7e-13 


29.5 


17.7 


19 


3.0 32 

— 


— 


— 
— 


— 
— 


1.9e-12 
6.3e-13 


30.5 
30 


379.4 
40.4 


322 
37 


- 0.15 

3*0.32 


— 


— 


— 


— 


2.7e-14 


28.7 


5.5 


6 


1^0.48 


— 


— 


— 


— 


2e-13 


29.5 


48.5 


45 


^0.33 
j|0.75 


— 


— 


— 


— 


4.5e-14 


28.9 


21.8 


27 










1.6e-13 


29.4 


48.8 


42 


if 0.22 


- 


- 


- 


- 


4.8e-14 


28.9 


52.8 


32 


S; 0.012 


— 


— 


— 


— 


5.4e-14 


29 


21.6 


16 


^ 0.16 










1.9e-14 


28.5 


0.6 


3 


^0.9 


- 


- 


- 


- 


2.3e-13 


29.6 


2 


6 


0.92 


- 


- 


— 


— 


2.6e-13 


29.7 


83.2 


75 


0.24 


— 


jj.tj 534| 

— 


— 
— 


— 
— 


3e-12 
l.le-14 


30.7 
28.3 


245.8 
2.8 


209 
8 


0.041 
0.95 




— 


— 


— 


3.9e-13 


29.8 


443.7 


474 


0.84 


, . _ „n so 

73 u ° u 

— 


52 20 52 26 


— 


— 

— 


1.9e-13 
2.8e-14 


29.5 
28.7 


127.7 
15.8 


100 
12 


0.032 
0.2 










O la 1 O 

z.je-lz 


5U.0 


1 co 1 n 


1 10 1 


U.UU /B 










l.le-13 


29.3 


121 


107 


0.17 










4e-14 


28.8 


0.6 


2 


0.76 










1.3e-13 


29.4 


4.4 


5 


0.5 










1.2e-14 


28.3 


6.2 


7 


0.52 










3.5e-14 


28.8 


5 


6 


0.54 



2.6; 
1.3 
1.6 
0.7 
4.5 
2.5 
1.7 



B 

,1:9 

AS 
0.6 

;3.2 

u 



13.4 
J 2.7 



1.6 



,46 

3-1 
■ 2.0 



4.6' 
4.5 
2.0 

!:?| 

3 6 7 - 

2 5 120 

4» J 

1 5 L6 
9 1 
19° 29 
1.5 1 - 6 

47^- 57 
"■^'0.42 
O92.0 

26" 

1 2 L5 
1 5 4 - 2 

c olO.O 
J.J 24 

27 140 
"■^'o.oi 

1 4 3 ' 



3.1 

2.6i 



q 7 7.4 

> 10 

5 9 64 

5-5 

-2 -24.4 
•'••'2.6 

2.2 
3 9° 
8.7 
9.1 



52 72' 2 : ^ 

J ,/ 52.68 

52 04 52 16 

51 47 5158 

J1.H/ 51 27 

53.36::;, 

52 13 5222 

J 52.05 

51.75g 247 
52.17= 2 - 79 
51 05 5169 
515? m 

JI,J1 SI .33 



Notes. 

Source 38 (Elias 29): Fit with variable abundance (Z/Z COf// > = 3.2 ~ Z Q ). 

Source 53 (Srl2A):P(^ 2 > xfy is low but the overall spectrum shape is well modeled with 3 thermal components. Abundance was a free parameter: Z = 0.3Zcoup 
Source 54: spectrum contamination from source nr. 53, see Sect. l3.2.2l for details. 

Source 61: the best fit model is the sum of two differently absorbed APEC models, the parameters in the table refer to the hot component. The soft component has parameters: N H 
kT = 0A7° f 3 , log E.M. = 50.47^ cnT 3 , absorbed flux: f x (0.3 - 10. keV band) = 1.67 • 10~ 15 erg cnT 2 s" 1 , unabsorbed flux: f x (0.3 - 10. keV band) = 1.99 • 10~ 15 erg cm 
luminosity: L x = 3.4 ■ 10 27 erg s -1 . 



I. Pillitteri et al.: Results from DROXO. , Online Material p 8 



Table A.3. Upper limits to count rate, X-ray flux and luminosity in the 0.3-10 keV fo r the ISOCAM YSOs dBontemps et all 1200 ll) undetected 
in the DROXO field of view. Count rates are scaled to MOS detector units as in Table |A~T1 (see Sect. l3.lt . Class III candidates from BoOl are 
indicated by "III?". Fluxes and luminosities are calculated by assuming a conversion factor from count rates of cf = 6.1 • 10~" and a distance of 
120 pc. For ISOCAM nr. 92 (WL 16) we have taken into account the large absorption toward this star (Ay ~ 30 mag, cf. Sect. 6.2) 



R.A. 


Dec. 


ISOCAM nr. 


Class 


Lim. rate 


Exposure time 


log F x 


log L x 


J2000 


J2000 






ct/ks 


ks 


1 9 

erg s cm - 


erg s _I 


16:26:30.9 


-24:31:07 


47 


III? 


0.798 


114.8 


4.87e-14 


28.92 


16:26:40.7 


-24:30:53 


55 


III? 


0.563 


186.3 


3.44e-14 


28.77 


16:26:41.6 


-24:40:15 


56 


11 


0.292 


350.6 


1.78e-14 


28.49 


16:26:42 


-24:33:24 


58 


III 


0.393 


267.3 


2.4e-14 


28.62 


16:26:42.1 


-24:31:02 


59 


11 


0.456 


225 


2.78e-14 


28.68 


16:26:52.1 


-24:30:39 


75 


11 


0.41 


254.7 


2.5e-14 


28.63 


16:26:53.6 


-24:32:36 


76 


II 


0.348 


305.8 


2.12e-14 


28.56 


16:26:56.7 


-24:28:38 


82 


III? 


0.514 


188.6 


3.14e-14 


28.73 


16:26:58.3 


-24:37:40 


85 


II 


0.328 


300.9 


2e-14 


28.54 


16:27:02.5 


-24:37:30 


92 


11 


0.222 


490.7 


2.31e-13 


29.60 


16:27:04.1 


-24:28:33 


95 


11 


1.163 


67.7 


7.09e-14 


29.09 


16:27:05.4 


-24:36:31 


99 


I 


0.218 


481.6 


1.33e-14 


28.36 


16:27:05.7 


-24:40:12 


100 


III? 


0.214 


470.4 


1.31e-14 


28.35 


16:27:07.9 


-24:40:27 


104 


III? 


0.192 


581.6 


1.17e-14 


28.31 


16:27:09.6 


-24:29:55 


109 


III? 


0.413 


239.6 


2.52e-14 


28.64 


16:27:10.3 


-24:33:22 


111 


III? 


0.308 


329 


1.88e-14 


28.51 


16:27:12.1 


-24:34:48 


115 


II 


0.265 


371.1 


1.62e-14 


28.45 


16:27:14.6 


-24:26:55 


118 


11 


0.803 


57.3 


4.9e-14 


28.93 


16:27:15.7 


-24:26:46 


120 


II 


1.953 


38.8 


1.19e-13 


29.31 


16:27:24.3 


-24:41:46 


136 


III? 


0.215 


651.6 


1.31e-14 


28.35 


16:27:24.8 


-24:41:03 


137 


I 


0.224 


651.6 


1.37e-14 


28.37 


16:27:26.2 


-24:42:45 


139 


II 


0.244 


340.4 


1.49e-14 


28.41 


16:27:36.3 


-24:28:34 


158 


III? 


0.672 


74.1 


4.1e-14 


28.85 


16:27:41.8 


-24:46:45 


170 


II 


0.255 


462.1 


1.56e-14 


28.43 


16:27:41.9 


-24:43:37 


171 


II 


0.222 


497 


1.35e-14 


28.37 


16:27:43.7 


-24:43:07 


173 


III? 


0.231 


441.8 


1.41e-14 


28.39 


16:27:45.9 


-24:37:60 


174 


III? 


0.215 


482.4 


1.31e-14 


28.35 


16:27:46 


-24:44:52 


175 


II 


0.236 


469 


1.44e-14 


28.39 


16:27:46.2 


-24:31:40 


176 


11 


0.5 


122.3 


3.05e-14 


28.72 


16:27:50 


-24:44:15 


179 


III? 


0.263 


418.4 


1.6e-14 


28.44 


16:27:50.3 


-24:39:01 


181 


III? 


0.542 


247.7 


3.31e-14 


28.76 


16:27:57.8 


-24:36:01 


189 


III? 


0.287 


350.3 


1.75e-14 


28.48 


16:28:05.5 


-24:33:55 


192 


III? 


0.367 


266.5 


2.24e-14 


28.59 


16:28:16.8 


-24:37:04 


196 


11 


0.963 


49.2 


5.88e-14 


29.01 


16:28:22.1 


-24:42:49 


197 


11 


0.848 


93.5 


5.17e-14 


28.95 



Table A.4. Spitzer IRAC fluxes of IR counterparts to DROXO X-ray sources from the catalog bv lEvans et al. (2003b). The a index is fitted from 
K to Spitzer MIPS1 band (2.2 - 24yum). The classification given in the last column is that given in lEvans et al.l d2003bl) . 



DROXO src. num. 


C2D id 


Name 


a 


IRAC Chi 


Ch2 


Ch3 


Ch4 


Mips 1 


Mips 2 


Classification 










3.6/jm 


4.5 fim 


5.8 yum 


8.0 pm 


24 urn 


70 yum 












mjy 


mjy 


mJy 


mjy 


mjy 


mjy 




1 


J162619.5-243727 




-1.35 


37.6 


28.8 


23.3 


26.3 


24.9 


_ 


YSOc: star+dust(IR4) 


2 


J162621.9-244440 


GY3 


-0.990 


28.8 


26.6 


23.8 


26.4 


33.6 





YSOc: star-t-dust(IRl) 


3 


J162623.7-244314 


DoAr25/GY17 


-1.12 


367. 


292. 


299. 


258. 


399. 


1100.0 


YSOc: star+dust(IR3) 


4 


J162627.5-244154 


GY33 


-0.510 


72.2 


59.0 


49.0 


50.9 


403. 


783. 


YSOc: star+dust(IR4) 


6 


J162635.2-244236 




-2.02 


0.123 


0.0967 


0.0794 


-0.0121 


_ 


_ 


two 


6 


J162635.3-244239 




0.540 


0.721 


1.17 


1.84 


2.82 


3.56 


_ 


Gale 


8 


J162644.2-243448 


WL12/GY111 


2.49 


239. 


744. 


1610.0 


2240.0 


— 


8120.0 


rising 


9 


J162644.3-244314 


GY112 


-2.59 


51.8 


34.8 


24.9 


14.6 


1.95 


_ 


star 


10 


J162644.4-244714 




-2.45 


21.0 


14.2 


10.0 


6.18 


1.35 


_ 


star 


11 


J162645. 1-245205 




0.310 


0.0496 


0.0676 


0.0703 


0.0676 


0.263 


_ 


two 


13 


J162647.0-244430 


GY122 


-2.43 


25.2 


17.4 


12.5 


7.26 


0.887 


_ 


star 


14 


J 162648. 1-244203 




-2.40 


24.0 


17.0 


12.3 


7.28 


0.875 


_ 


star 


15 


J162648.5-242839 


WL2/GY128 


0.02 


79.2 


111. 


128. 


148. 


446. 


1100.0 


YSOc: red 


17 


J162649.0-243825 


WL18/GY129 


-0.930 


101. 


135. 


103. 


86.7 


122. 


_ 


YSOc 


18 


J162652.6-244332 




-2.54 


0.256 


0.178 


0.123 


0.0706 


0.455 


_ 


star 


18 


J162652.9-244324 




-1.04 


0.0511 


0.0506 


0.0109 


0.0422 


_ 


_ 


two 


21 


J162656.5-244960 




0.4 


0.0628 


0.111 


0.183 


0.206 


0.304 


_ 


cup-down 


22 


J162657.3-243539 


WL21/GY164 


-0.660 


7.70 


12.3 


11.9 


10.7 


12.1 


_ 


YSOc 


23 


J162657.5-244606 




-2.67 


32.8 


20.7 


14.3 


8.49 


_ 


_ 


star 


24 


J162657.6-244313 




-0.790 


0.0505 


0.0530 


0.0716 


0.0152 


_ 


_ 


two 


25 


J162658.5-244537 


SR24 


-0.920 


990. 


-30.6 


351. 


1040.0 


1320.0 


_ 


star+dust(IR4) 


25 


J162658.4-244532 




-0.450 


1390.0 


1440.0 


1740.0 


1950.0 


2230.0 


6380.0 


YSOc: star+dust(IRl) 


27 


J162659.2-243459 


WL22/GY174 


2.36 


143. 


487. 


564. 


3230.0 


4540.0 


_ 


rising 


28 


J162659.0-243557 


WL14/GY172 


-1.01 


14.2 


13.7 


11.6 


10.4 


20.0 


_ 


YSOc: star+dust(IR4) 


28 


J162659.0-243557 


WL14/GY172 


-1.01 


14.2 


13.7 


11.6 


10.4 


20.0 


_ 


YSOc: star+dust(IR4) 


30 


J162704.5-244260 


GY193 


-2.46 


48.8 


35.1 


25.3 


15.3 


1.79 


_ 


star 


31 


J162704.5-244214 


GY194 


-2.45 


48.6 


36.3 


26.2 


15.4 


1.79 


_ 


star 


32 


J162705.2-244113 




1.07 


0.0330 


0.0946 


0.160 


0.177 


0.279 


_ 


cup-down 


33 


J162706.6-244149 


GY204 


-0.870 


29.4 


24.3 


21.6 


22.3 


57.9 


_ 


YSOc: star+dust(IR2) 


34 


J162706.8-243815 


WL17/GY205 


0.610 


240. 


416. 


553. 


695. 


2790.0 


6070.0 


YSOc: red 


35 


J162709. 1-243408 


WL10/GY211 


-0.910 


259. 


310. 


272. 


222. 


339. 


784. 


YSOc: star+dust(IR2) 


36 


J162709.3-244320 




-2.13 


21.8 


17.1 


13.1 


7.79 


0.858 


_ 


star 


37 


J162709.3-244022 


GY213 


-0.09 


44.1 


53.7 


65.1 


92.9 


202. 


_ 


YSOc: star+dust(IR3) 


39 


J16271 1.2-244047 


GY224 


-0.05 


203. 


305. 


358. 


367. 


907. 


908. 


YSOc 


40 


J16271 1.7-243832 


WL19/GY227 


-0.430 


215. 


354. 


406. 


328. 


223. 


— 


YSOc 


42 


J162713. 8-244332 


GY235 


-0.220 


56.2 


60.7 


70.8 


116. 


309. 


577. 


YSOc: star+dust(IR2) 


43 


J162715. 1-245139 




-1.06 


82.5 


72.8 


58.6 


52.3 


139. 


_ 


YSOc: star+dust(IR2) 


43 


J162714.5-245133 




-2.53 


26.4 


17.3 


11.9 


7.24 


1.28 


_ 


star 


44 


J1627 15.5-243054 


IRS35/GY238 


0.010 


19.0 


31.4 


41.1 


50.7 


75.4 


- 


YSOc 


46 


J1627 15.7-243843 


WL20/GY240 


-0.430 


82.1 


13.0 


130. 


59.2 


993. 




red 


46 


J162715.9-243843 


WL20/GY240 


-0.7 


127. 


143. 


140. 


99.0 


816. 


15700. 


star+dust(IR2) 


46 


J1627 15.9-243843 


WL20/GY240 


-0.7 


127. 


143. 


140. 


99.0 


816. 


15700. 


star+dust(IR2) 


47 


J1627 16.4-243 114 




0.160 


19.7 


52.7 


66.5 


50.7 


129. 




YSOc 


49 


J162718.2-242853 


WL5/GY246 


-0.520 


209. 


297. 


298. 


163. 






cup-down 


50 


J 1627 18. 4-245454 




-1.24 


-999. 


51.2 


-999. 


38.2 


39.0 




star+dust(IR2) 


51 


J162718.4-243915 


GY245 


0.180 


43.8 


72.0 


91.3 


107. 


251. 




YSOc 



Table A.4. continued. 



DROXO src. num. 


C2D id 


Name 


a 


IRAC Chi 


Ch2 


Ch3 


Ch4 


Mips 1 


Mips 2 


Classification 










3.6 yum 


4.5 jum 


5.8 /im 


8.0 pm 


24 pm 


70 fim 














mjy 


mJy 


mJy 


mJy 


mJy 


mJy 






52 


J162718.5-242906 


WL4/GY247 


-0.590 


173. 


195. 


232. 


252. 


269. 


- 


YSOc: 


star+dust(IRl) 


52 


J162719.2-242844 


WL3/GY249 


-0.03 


98.6 


163. 


208. 


204. 


325. 


- 


YSOc 




53 


J162719.5-244140 


SR12/GY250 


-2.56 


140. 


94.7 


68.3 


39.9 


5.24 


- 


star 




54 


J162721.5-244143 


IRS42/GY252 


-0.03 


1060.0 


1630.0 


2100.0 


2980.0 


3450.0 


2940.0 


YSOc 




55 


J162721. 8-244336 


GY253 


-1.89 


64.7 


60.4 


53.3 


30.9 


3.18 


- 


star 




56 


J162721. 8-242953 


WL6/GY254 


0.720 


467. 


925. 


1440.0 


1730.0 


4360.0 


5110.0 


YSOc 




57 


J162723.0-244807 




-1.04 


101. 


89.1 


84.8 


104. 


91.7 


- 


YSOc: 


star+dust(IR2) 


59 


J162724.6-242935 




-1.45 


13.0 


12.4 


10.9 


6.42 


4.26 


- 


YSOc: 


star+dust(MPl) 


60 


J162726.5-243923 


GY262 


-0.460 


182. 


237. 


222. 


210. 


257. 


- 


rising 




61 


J 162727. 1-2432 17 




-1.54 


11.6 


10.7 


9.18 


5.37 


1.04 


- 


star 




62 


J162726.9-244051 


IRS43/GY265 


1.17 


629. 


1240.0 


1790.0 


2190.0 


- 


34400. 


rising 




63 


J162727.4-243117 


VSSG25/GY267 


-0.940 


130. 


115. 


101. 


97.6 


196. 


- 


YSOc: 


star+dust(IR4) 


64 


J162728.0-243933 


IRS44/GY269 


2.29 


731. 


1830.0 


2940.0 


2320.0 


- 


34700. 


rising 




65 


J 162728. 4-242721 


IRS45/GY273 


-0.03 


187. 


272. 


382. 


481. 


712. 


8790.0 


YSOc 




66 


J162728.7-245432 




-2.40 


-999. 


5.27 


-999. 


2.24 


0.923 


- 


star 




67 


J162729.4-243916 


IRS46/GY274 


0.180 


172. 


271. 


402. 


411. 


639. 


- 


rising 




68 


J162730.2-242743 


IRS47/GY279 


-0.120 


740. 


1190.0 


1580.0 


2040.0 


1720.0 


7040.0 


YSOc 




69 


J162729.9-243336 


GY278 


-1.38 


22.1 


21.6 


19.6 


11.2 


0.865 


- 


YSOc 




70 


J162730.5-245211 




-2.41 


0.0755 


0.0541 


0.0266 


0.103 


0.0906 


- 


two 




70 


J162730.5-245208 




0.08 


0.256 


0.324 


0.404 


0.637 


1.74 


- 


Gale 




71 


J162730.8-244727 


B 162730-244726 


-2.40 


80.9 


60.9 


47.7 


26.7 


2.22 


- 


star 




72 


J16273 1.0-243403 


GY283 


-2.34 


35.6 


26.8 


19.8 


11.7 


1.51 


- 


star 




73 


J16273 1.3-243 110 




0.720 


0.0933 


0.163 


0.187 


0.290 


1.18 


- 


red 




74 


J162732.7-2445 




-1.48 


6.26 


6.05 


4.98 


3.10 


1.97 


- 


YSOc: 


star+dust(MPl) 


75 


J162732.6-243323 


GY289 


-0.920 


45.6 


41.9 


36.8 


32.8 


56.5 


- 


YSOc: 


star+dust(IR4) 


76 


J162732.8-243235 


GY291 


-0.460 


54.2 


69.5 


84.8 


108. 


94.5 


- 


YSOc: 


star+dust(IR2) 


77 


J162733.1-244115 


GY292 


-0.880 


571. 


494. 


473. 


558. 


721. 


- 


YSOc: 


star+dust(IR4) 


78 


J162735.3-243833 


GY295 


-2.47 


54.8 


39.3 


27.9 


17.0 


1.69 


- 


star 




79 


J162735.7-244532 


GY296 


-2.28 


22.0 


15.6 


11.3 


6.92 


0.507 


- 


star 




80 


J162737.2-244238 


GY301 


0.130 


91.6 


142. 


162. 


179. 


561. 


845. 


YSOc: 


red 


81 


J162738. 1-243043 




-1.89 


72.0 


53.8 


43.3 


13.5 


- 


- 


star 




82 


J162738.3-243658 


IRS49/GY308 


-0.730 


375. 


335. 


331. 


397. 


688. 


924. 


YSOc: 


star+dust(IR3) 


83 


J162738.6-243839 


GY310 


-0.490 


23.7 


22.3 


21.2 


29.2 


87.6 


- 


YSOc: 


star+dust(IR2) 


84 


J162738.9-244020 


GY312 


0.640 


15.8 


26.7 


38.2 


53.2 


461. 


950. 


YSOc: 


star+dust(IRl) 


85 


J162739.0-244721 




0.4 


0.170 


0.265 


0.430 


0.602 


1.63 


- 


Gale 




86 


J162739.4-243915 


GY314 


-0.370 


282. 


284. 


333. 


427. 


1140.0 


1310.0 


YSOc: 


star+dust(IR2) 


87 


J162739.8-244315 


IRS51/GY315 


-0.150 


752. 


916. 


1000.0 


1070.0 


2730.0 


3260.0 


YSOc 




89 


J162741.5-243538 


GY322 


-1.98 


38.5 


33.2 


24.7 


14.7 


5.77 


- 


YSOc: 


star+dust(MPl) 


90 


J162742.7-243851 


GY326 


-1.28 


31.4 


25.4 


21.4 


20.1 


24.5 


- 


YSOc: 


star+dust(IR4) 


91 


J162743.6-243127 




1.72 


0.0248 


0.0560 


0.123 


0.175 


- 


- 


cup-down 


92 


J162744.4-244856 




-1.33 


0.0636 


0.0589 


0.0455 


0.0251 


0.136 


- 


two 




93 


J162747. 1-244535 


GY352 


-0.720 


41.1 


47.8 


52.5 


49.9 


57.1 




YSOc: star+dust(IR2) 


94 


J162749.1-243906 




0.510 


0.187 


0.357 


0.614 


0.908 


1.70 




Gale 




97 


J162751.8-243145 


IRS54/GY378 


0.03 


525. 


712. 


931. 


1010.0 


3560.0 


6500.0 


YSOc: red 


98 


J162751.9-244630 


GY377 


-2.26 


40.3 


31.6 


25.4 


14.6 


1.44 




star 




99 


J162752.1-244050 


IRS55/GY380 


-2.61 


203. 


145. 


106. 


62.7 


7.40 




star 




1 


J162755.6-244451 


GY398 


-2.36 


34.2 


26.0 


18.3 


10.7 


1.28 




star 





Table A.4. continued. 



DROXO src. num. 


C2D id 


Name 


a 


IRAC Chi 


Ch2 


Ch3 


Ch4 


Mips 1 


Mips 2 


Classification 










3.6 yum 


4.5 jum 


5.8 /im 


8.0 pm 


24 pm 


70 fim 












mjy 


mJy 


mJy 


mJy 


mJy 


mJy 




101 


J162757. 8-244002 


GY410 


-2.42 


54.6 


40.3 


30.0 


17.9 


1.69 


- 


star 


102 


J162759.9-244819 




-2.56 


72.8 


52.3 


36.7 


22.1 


2.16 


- 


star 


103 


J162760.0-244239 




0.660 


0.0625 


0.102 


0.159 


0.209 


0.415 


- 


cup-down 


104 


J162802.0-244953 




-1.14 


0.103 


0.101 


0.0904 


0.105 


0.441 


- 


star 


104 


J162802.1-244956 




-2.60 


0.113 


0.0776 


0.0361 


0.0274 


0.356 


- 


two 


105 


J162804.6-243456 


GY463 


-1.85 


52.7 


48.2 


41.0 


24.7 


3. 


- 


star 


106 


J162804.8-243710 




-2.13 


43.2 


33.4 


25.6 


14.3 


4.72 




YSOc: star+dust(MPl) 


107 


J162804.9-245014 




-2.73 


0.0902 


0.06 


0.00124 


0.0493 


0.524 




two 


107 


J162805.3-245019 




-2.34 


0.0542 


0.0395 


0.00871 


0.000535 


0.303 




two 


108 


J162806.8-244828 




0. 


0.405 


0.446 


0.856 


0.989 


1.93 




Gale 


109 


J162808. 1-244512 




-2.22 


6.38 


4.79 


3.41 


2.18 


0.303 




star 


111 


J162810.2-243928 




-2.04 


0.267 


0.209 


0.0903 


0.0890 






two 
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fective temperatures and bolometric luminosities of ISOCAM objects in DROXO (see Sect. 13.31 1. 



ISO 


DROXO Src. 


Class 


log T eff (K) 


log L/L bo i 


log M/MO 


1 




II 


3.66 


0.3 


-0.18 


2 




11 


3.59 


-0.05 


-0.41 


3 




n 


3.61 


0.03 


-0.35 


5 




in 


3.65 


0.26 


-0.2 


6 




ii 


3.64 


0.21 


-0.23 


9 




ii 


3.47 


-0.92 


-0.91 


11 




in 


3.57 


-0.13 


-0.49 


12 




ii 


3.48 


-0.86 


-0.88 


13 




ii 


3.58 


-0.1 


-0.45 


14 




in 


3.63 


0.15 


-0.26 


16 




in 


3.68 


0.61 


0.04 


17 




ii 


3.68 


0.65 


0.07 


18 




in 


3.64 


0.2 


-0.23 


19 




ii 


3.68 


0.67 


0.08 


20 




ii 


3.63 


0.13 


-0.28 


23 




ii 


3.43 


-1.6 


-1.32 


24 




ii 


3.67 


0.39 


-0.12 


26 




ii 


3.54 


-0.31 


-0.62 


28 




in 


3.67 


0.48 


-0.06 


30 




ii 


3.44 


-1.2 


-1.07 


32 


2 




3.43 


-1.47 


-1.24 


33 






3.38 


-2.93 


-2.26 


35 




ii 


3.46 


-0.96 


-0.93 


36 




ii 


3.69 


0.81 


0.17 


37 






3.51 


-0.6 


-0.77 


38 


3 




3.63 


0.15 


-0.27 


39 




ii 


3.76 


1.51 


0.47 


40 




ii 


3.68 


0.65 


0.06 


41 






3.53 


-0.43 


-0.69 


43 


4 


ii 


3.63 


0.16 


-0.26 


44 




in 


3.6 





-0.36 


46 




ii 


3.61 


0.04 


-0.34 


47 




in 


3.55 


-0.23 


-0.57 


51 




ii 


3.52 


-0.49 


-0.72 


52 




ii 


3.51 


-0.6 


-0.77 


53 




ii 


3.51 


-0.62 


-0.78 


55 




in 


3.43 


-1.57 


-1.3 


56 




ii 


3.53 


-0.39 


-0.66 


58 




in 


3.68 


0.7 


0.1 


59 




ii 


3.55 


-0.23 


-0.57 


59 




ii 


3.57 


-0.13 


-0.48 


62 






3.68 


0.66 


0.07 


63 




ii 


3.51 


-0.55 


-0.75 


64 




in 


3.64 


0.19 


-0.24 


66 


9 


in 


3.52 


-0.48 


-0.72 


67 




ii 


3.66 


0.34 


-0.15 


68 




ii 


3.67 


0.43 


-0.09 


69 


13 


in 


3.49 


-0.75 


-0.83 


70 


15 


ii 


3.38 


-3.01 


-2.41 


72 


17 


ii 


3.5 


-0.67 


-0.8 


73 




in 


3.68 


0.68 


0.08 


74 




in 


3.58 


-0.09 


-0.45 


75 




ii 


3.49 


-0.71 


-0.82 


76 




ii 


3.52 


-0.5 


-0.72 


78 




ii 


3.66 


0.33 


-0.16 


79 




ii 


3.44 


-1.27 


-1.11 


82 




in 


3.43 


-1.73 


-1.4 


83 




ii 


3.62 


0.11 


-0.29 


84 


22 


11 


3.46 


-1.02 


-0.96 


84 






3.46 


-1.01 


-0.96 


86 






3.53 


-0.42 


-0.68 


87 






3.47 


-0.9 


-0.9 


88 


25 




3.67 


0.5 


-0.04 


88 






3.67 


0.53 


-0.02 


89 


28 




3.5 


-0.68 


-0.8 


89 






3.5 


-0.66 


-0.8 


91 




in 


3.69 


0.75 


0.13 
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Appendix B: Atlas of DROXO sources. 
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We sh ow here an example of the atlas that we have produced for each DROXO source. It is reported the number of the source, the literature 
name as in Bontemps et al. (2001), the EPIC image, the J-band 2MASS image, the EPIC spectrum, and the light curve. The atlas is available only 
online at the following web address: http:www.astropa.unipa.it/ pilli/atlas_droxo_sources.pdf. 
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